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ABSTRACT 



This report details the development and provides the 
documentation for custom computer software applications that 
evaluate antenna parameters. The applications are written in 
Mathcad 3.1 for the following antenna types: linear, planar, 
and circular arrays; folded dipole; caged dipole; parabolic 
reflectors with helical and spiral feeds. Inputs to the 
applications are limited to the antenna's physical dimensions. 
In some cases, ground parameters are required. 

The chapters are structured to provide the user with the 
necessary information needed to use and interpret the software 
for each antenna type. The software applications are provided 
as appendices and give examples of each antenna type. 

Outputs of the applications provide various numerical and 
performance predictions, as well as far-field radiation 
patterns. The results computed are consistent with 
predictions provided in applicable publications, as well as 
those calculated by numerical antenna analysis programs such 
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as ELNEC. 
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I. BACKGROUND AND PURPOSE 



This thesis details the development of user friendly 
mathematical applications capable of computing the radiation 
pattern and other pertinent antenna parameters of an antenna 
or antenna system based on available information, primarily 
dimensional information obtained from photographs. The 
applications developed are compatible with any IBM personnel 
computer using MS-DOS version 3.2 or higher and a math 
coprocessor. The applications are written in Mathcad version 
3.1 engineering software and can be run using either Mathcad 
3 . 1 or Mathcad version 4.0. These applications are intended 
for use by the Naval Maritime Intelligence Center 
(NAVMARINTCEN) . 

As already mentioned, user inputs are limited to antenna 
dimensions based on physical measurements primarily obtained 
from photographs, although other source data can be used. The 
applications are flexible to the extent that other parameters 
may be estimated and used as inputs to increase the accuracy 
of the computations to gain better insight into the antenna's 
performance. The Mathcad applications provide various 
performance predictions as well as a graphical representation 
of the antenna's far-field radiation pattern. The necessary 
background information needed to interpret the application's 
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formulas and displays are provided in the corresponding thesis 
chapters . 

Dietrich [Ref. 1] and Gerry [Ref. 2] have completed the 
first and second reports of this project, respectively. This 
thesis is the last in a series of three reports intended to 
fulfill the NAVMARINTCEN statement of work. 



2 



II. INTRODUCTION 



The increasingly sophisticated design by foreign 
countries of antennas warrants the development of improved 
analytical methods to provide timely and accurate technical 
assessment of these antennas. Through the process of reverse 
engineering, foreign communication, navigation, 
Identification, Friend or Foe (IFF) , and radar antennas are 
evaluated to determine their capabilities, limitations, and 
vulnerabilities. Having limited information available 
concerning the parameters of an antenna, usually restricted to 
photographs that provide only the type and physical dimensions 
of the antenna, current methods of antenna evaluation are slow 
and tedious. 

With the availability of high-speed personal computers 
matched with sophisticated off-the-shelf engineering software, 
the assessment of antennas using only photographic information 
may now be obtained rapidly. The objective of this report is 
to provide NAVMARINTCEN with user friendly Mathcad software to 
achieve this type of performance in their evaluation of 
antenna systems. 

As with any computer analysis of antennas, the accuracy 
of the results depend upon the sophistication of software 
used, the complexity of calculations, the reliability of input 
data obtained from photographs, and the accuracy of any 
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required input data that must be estimated. This report and 
the Mathcad applications are written with existing engineering 
equations for the type of antenna analyzed and inputs are 
limited to physical dimensions. Additional information such 
as ground effects and characteristic impedance may be included 
to increase the accuracy of the analysis. To this end, this 
report should provide an initial understanding as to the 
capabilities of the antenna system in question. 

The chapters of this report document the equations and 
assumptions used for each specific antenna type and is written 
as a comprehensive reference for the software. Copies of each 
application are included as appendices to provide the user 
with a printed illustration of the software. 
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III. THE UNIFORM LINEAR ARRAY ANTENNA 



Arrays are highly directional antenna system with narrow, 
steerable beams and low side lobes used for long distance 
communication and radar systems. The linear array is one 
geometrical configuration that consists of radiating elements 
lying along a straight line. As a result of the dependance of 
inter-element spacing (d) to the wavelength (X) , linear arrays 
are inherently narrowband antenna systems. 

In this Mathcad application uniformly spaced, equally 
excited, isotropic point sources aligned on the z-axis are 
used as the radiating elements. In using isotropic point 
sources, polarization of the array is not calculated and an 
assumption is made that the antenna input resistance (R in ) is 
equal to the radiation resistance (R r ) . This assumption 
results in calculations for gain (G) that are idealized 
values. The following inputs are needed to implement the 
software : 

N = number of isotropic radiating elements 
d = inter-element spacing between adjacent elements 
0 O = direction of main lobe 
f - frequency of interest 
I c = antenna feed current 
Z 0 = characteristic feed impedance 
The first two inputs are obtained from visual data such as 
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photographs, and the * indicates input parameters that are 
either known or are estimated. The linear array geometry is 
illustrated in Figure 3.1. 




FIGURE 3.1 Geometry of N-element linear array of 
isotropic point sources 
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The total electric field of a linear array in the far- 
field, neglecting mutual coupling between adjacent elements, 
is determined by the vector addition of the fields of 
individual radiating elements. It can be shown that the total 
field is equivalent to the product of the field of a single 
element, selected at a reference point, and a factor which is 
referred to as the array factor (AF) : 

E( total) = [E (element)] x [AF] (3.1) 

In (3.1), E is the vector electric field intensity, while AF 
is a scalar quantity. This result is referred to as pattern 
multiplication for arrays of identical elements. The array 
factor is a function of the geometry of the array and the 
element excitation amplitude and phase. For uniform 
excitation, by varying the inter-element spacing (d) and/or 
the phase (p) between the elements, the characteristics of the 
array factor, and consequently the total electric field, of 
the array can be controlled. [Ref. 3: pp. 204-207] 

This Mathcad application is written for an N-element 
linear array with equally spaced, uniformly excited isotropic 
radiating elements. In addition, element phase is assumed to 
vary linearly along the length of the array. The radiation 
pattern of a specific element type is neglected in (3.1) since 
in normal usage it will have little effect when the array 
consists of a large number of elements. However, the user of 
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this software should be able to use information regarding the 
type of element and the orientation of the elements to narrow 
the range of possible 0 o 's. The normalized array factor (AF n ) 



is [Ref. 3: pp. 212-214]: 




( dimensionless ) (3.2) 



where 



\J/(0) = kdcos(Q) + P ( radians ) 



(3.3) 



P = -kdcos(Q 0 ) ( radians ) 



(3.4) 



In (3.4), 0 O is the direction of the maximum array factor 
value (i.e., for a broadside array 0 O = n/2, for an end-fire 
array with the main beam directed at 180°, 0 O = n , and for a 
phased array with the main beam directed to ±60°, 0 O = 7t/3), 
P is the phase excitation difference between the elements, k 
is the wavenumber (2n/X) , and d is the inter-element spacing. 

With the linear array aligned along the z-axis, the array 
factor is a function of 0 but not of 0. The visible region of 
the array in terms of 0, i(i(0), and the inter-element spacing 
(d/A) is: 



0 < 0 < 7i ( radians ) 



(3.5) 
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$-kd < (0) < $ + kd ( radians ) 



(3.6) 



d/A determines how much of the array factor appears in the 
visible region as defined by (3.5) and (3.6). Although d/A 
provides guidance in determining the frequency to select for 
this Mathcad application, it should be noted that most arrays 
are designed with the inter-element spacing less than one 
wavelength, usually close to a hal f-wavelength . Thus, linear 
arrays are characteristic narrowband antenna systems. [Ref. 
4: p. 123] 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 

r k 1.6 A ( meters ) (3.7) 

i k S (Nd) ( meters ) (3.8) 



r * 2 W ) 2 



( meters ) 



(3-9) 



The array length, Nd, in (3.8) and (3.9) is the maximum 
dimension of the antenna (the array length is assumed to 
include a distance of d/2 beyond each end element [Ref. 5: 
p . 55 ] ) . The minimum distance to the far-field is found by 

selecting the largest value of (3.7), (3.8), and (3.9). 

Radiation intensity (U(6)) is related to the power 
radiated in a given direction and is independent of the 
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distance to the observation point [Ref. 3: p. 27; Ref. 4: p. 
33]. For linear arrays U(0) is given by [Ref. 3: pp. 229- 
233] : 

[7(0) = ( AF n ) 2 (W / solid angle) (3.10) 

Radiated power (P rad ) is the total power radiated by the 
antenna and is obtained by integrating U(0) over a surface 
surrounding the antenna [Ref. 4: p. 33]. Since the geometry 
of the linear array is aligned on the z-axis and the array 
factor is independent of 0, P rad is thus [Ref. 3: pp. 229-233 ]: 

P rad = 2n f n [/(d) sin ( 0 ) c?0 (N) (3.11) 

J 0 

Directivity (D 0 ) is the maximum value of directive gain 
where directive gain is the ratio of U(0) in a specific 
direction to the average radiation intensity [Ref. 4: pp. 34- 
36]. Hence, D 0 is given by [Ref. 3: pp. 229-233]: 

4 C7 m . v 

D o = — ( dimensionless ) (3.12) 

° i ad 

In (3.12), U max equals unity and occurs at 0 O for the uniform 
linear array. 

Effective isotropic radiated power (EIRP) is a commonly 
used communications term that is defined as the product of 
antenna gain and total power accepted by the antenna from the 
transmitter. EIRP is determined as [Ref. 4: p. 62]: 
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EIRP = P 



(fv T ) 



(3.13) 



Assuming that the power delivered to the array consisting 
of N-isotropic point sources is equivalent to the radiated 
power, we can estimate the radiation resistance (R r ) of the 
antenna as [Ref. 3: p. 55]: 



Using the same logic, we obtain the input resistance (R in ) of 
the array made up of isotropic point sources as equivalent to 
the radiation resistance: 



The exact gain (G) of the array cannot be calculated with 
the data assumed, but an idealized value can be computed based 
on the assumption made to determine R in and the assumptions 
made to calculate total antenna efficiency (e t ). The total 
antenna efficiency is given by [Ref. 3: pp. 44-45]: 




(3.14) 



Rin = Rr 



(3.15) 



e c = t r t cd = 1- r| 2 ( dimensionless ) 



(3.16) 



where 



t cd = 1 ( dimensionless ) 



(3.17) 



e 



I 



= i - 1 r | 2 



( dimensi onl ess) 



(3.18) 
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( dimensionless ) 



(3.19) 




O 



o 



In (3.17), e cd is the conduction and dielectric efficiency 
which is unity for an ideal lossless antenna. In (3.19), Z Q 
is the characteristic impedance of the transmission line and 
r is voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(3.18) and the ideal gain for a linear array is [Ref. 3: pp. 
43-44 ] : 



As a result of constructing the linear array with 
isotropic point sources, the polarization loss factor (PLF) 
and antenna polarization are not calculated. Therefore, the 
maximum effective aperture (A em ) of the array when used as a 
receiver is calculated assuming negligible polarization 
mismatches. The A era is then [Ref. 3: p. 63]: 



G = t t D 0 ( dimensionless ) 



(3.20) 



G(dB) = io log 10 ( t c D 0 ) (dB) 



(3.21) 



em 




(3.22) 



where 



PLF = 1 (dimensionless) 



(3.23) 
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The assumptions used in finding R r and A em result in an 
idealized value in calculating the effective height (h em ) . The 



effective height 



is [Ref. 



Km = 



6 : p . 



KKm 



42] : ] 

(meters) 



(3.24) 



where r) 0 is the characteristic impedance of free space (120 tt) . 

Two displays are produced in this Mathcad application. 
The first display is the rectangular representation of the 
magnitude of the array factor from -n to n . The second 
display is the polar representation of the magnitude of the 
array factor. 

The Hansen-Woodyard end-fire array was not addressed in 
this application nor were arrays of uniform spacing but 
nonuniform amplitude (i.e., Dolph-Tschebyschef f array). In 
addition, nonuniformly spaced, uniformly excited arrays were 
not addressed by this application. 

The results obtained with this application are identical 
to those found in [Ref. 3: pp. 216-240] with respect to array 
factor patterns and directivity. Table 3.1 is a comparison of 
data from [Ref. 3] for a 10-element equally excited, uniformly 
spaced, end-fire linear array while Figure 3.2 is a comparison 
of the array factor pattern of the same array. Of note, the 
array factor patterns from this application and the patterns 
from [Ref. 3] are identical as expected. 
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Table 3.1 End-fire Array Data Comparison 
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FIGURE 3.2 Array factor pattern of 10-element equally 
excited, uniformly spaced end-fire linear array with d = 
X/4 m, X 0 = 0°, and /3 = -kd 
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Table 3.2 is a comparison of data from [Ref. 3] for a 10- 
element equally excited, uniformly spaced, broadside linear 
array, and Figure 3.3 is a comparison of the array factor 
pattern of the same array. The array factor patterns from 
this application and the patterns from [Ref. 3] are identical. 



Table 3.2 Broadside Array Data Comparison 



ANTENNA 


REFERENCE 


CALCULATED 


DIMENSIONS 


DIRECTIVITY (D 0 ) 


DIRECTIVITY (D 0 ) 


N = 10 






d = A/4 m 


6.99 dB 


7.14 dB 


0 O = tt/2 






O 

II 
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FIGURE 3.3 Array factor pattern of 10-element equally 
excited, uniformly spaced broadside linear array with d = 
A/4 m, 0 O = 7r/2, and (3 = 0 
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IV. THE PLANAR ARRAY ANTENNA 



Planar arrays, frequently used long distance 
communications and radar systems, exhibit characteristics 
analogous to those of linear arrays but have additional 
capability to control and shape the radiation pattern. Planar 
arrays are more versatile than linear arrays, providing 
symmetrical patterns with lower side lobes, and can be used to 
scan the main beam of the antenna array toward any point in 
space [Ref. 3: p. 261]. As with linear arrays, planar arrays 
are narrowband antenna systems due to the dependance of the 
inter-element spacing (d x and d y ) on the wavelength (X) . The 
following inputs are required analyze a planar array: 

M = number of isotropic radiating elements in the 
x-direction 

N = number of isotropic radiating elements in the 
y-direction 

d x = inter-element spacing between adjacent elements in 
the x-direction 

d y = inter-element spacing between adjacent elements in 
the y-direction 
f = frequency of interest * 

0 O = direction of main lobe at 0 = 0 O * 

<P Q = direction of main lobe at </> = <p 0 " 

I Q = antenna feed current * 
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Z 0 = characteristic feed impedance 
The first four inputs are data obtained primarily from visual 
data such as photographs, and the * indicates input parameters 
that are either known or estimated. The planar array geometry 
is illustrated in Figure 4.1. 



z 



A 



e 




/ 



v 



X 




FIGURE 4.1 Geometry of M x N element planar array 
of isotropic point sources 
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This Mathcad application is written for an M x N planar 
array with uniformly spaced, equally excited, identical 
radiating elements positioned on the x-y plane. This 
application may also be used to approximate the array factor 
and field patterns for arrays that have elliptical or circular 
geometries by assuming an ellipse is equivalent to a rectangle 
or a circle is equivalent to a square, respectively [Ref. 6: 
p . 187 ] . Since isotropic point sources are assumed as array 
elements, the polarization of the array is not calculated, and 
an assumption is made that the antenna input resistance (R in ) 
is equivalent to the radiation resistance (R r ) . This 
assumption results in values for gain (G) that are idealized. 
As with the linear array, the radiation pattern of a specific 
element type is neglected in the analysis of the planar array 
since in normal usage it will have little effect for a planar 
array with a large number of elements. However, as with the 
linear array, the user of this software should be able to use 
information regarding the type of element and the orientation 
of the elements to narrow the range of the possible 0 o ,s and 
<Z> 0 's. 

The normalized array factor (AFJ is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [Ref. 3: pp. 
260-263 ] : 
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AF n (0,<|>) 



(4.1) 



sinj 






sinj 


[f V 0 .*)) 


M sin 


' tx( e '4>) 
2 




N sin 


' ty(0'4>) ) 

2 JJ 



where 

t x ( 0, 4> ) = kd x sin(0)cos(<{>) + p x ( radians) (4.2) 

4f y (0,<t>) = kd y sin(0) sin(<t>) + P y ( radians ) (4.3) 

P x = -kd x sin ( 0 O ) cos ( <t> 0 ) ( radians ) (4.4) 

p y = -Jcd y sin ( 0 O ) sin ( (J> 0 ) ( radians ) (4.5) 

Equation (4.1) is the normalized array factor that is a 
function of both 0 and <p where (4.2) and (4.3) are the arra\ 
factor phase shifts in the x and y directions, respectively. 
(4.4) and (4.5) are the progressive phase shift between 
adjacent elements in the x and y directions, respectively, 
with the main beam directed along 0 = 0 C and <p = <p 0 . In (4.2) 
through (4.5) , k is the wavenumber (2u/X) . 

When the inter-element spacing between the elements (d x 
and d y ) equals or is greater than the wavelength (A), multiple 
maxima of equal magnitude are formed. To avoid creating 
multiple maxima, the inter-element spacing must be less than 
a wavelength at the frequency of operation [Ref. 3: p. 263]. 
With the dependance of inter-element spacing on the 
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wavelength, planar arrays are characteristic narrowband 
antenna systems. 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 

r 2 : 1.6/1 (meters) (4.6) 



r 2 5 yj (Md x ) 2 + (Ndyp ( meters ) 



r 2 



2((Md x ) 2 ^ ( Nd ^ ) 2 ) 



(meters) 



(4.7) 



(4.8) 



The length of an array is assumed to include a distance d/2 
beyond each end element [Ref. 5: p. 55]. Therefore, the 
maximum dimension of the planar array is: 

^ ( Md x ) 2 + ( Wd y ) 2 (meters) (4.9) 



The minimum distance to the far-field is found selecting the 
largest value of (4.6), (4.7), and (4.8). 

Directivity (D 0 ) for planar arrays is determined 
approximately by [Ref. 7: p. 78]: 



D n - 4 7t 

^ 2n r n/2 



AF n AFn Lax 



f‘ (AF„AF' ) sin ( 0 ) d0 d(f> 



(W) 



(4.10) 



The maximum value of the numerator in (4.10) is unity at 0 = 
0 O and (p = <p 0 since the array factor is normalized. The 
denominator is the radiated power (P rad ) of the antenna where 
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P rad is obtained by integrating over all angles around the 
antenna above the x-y plane (0 < 0 < n / 2) [Ref. 3: p. 28; Ref. 
7: p. 78]. This assumes the array is always pointed toward 
one side or the other of the x-y plane. 

Effective isotropic radiated power (EIRP) for planar 
arrays is defined as the product of antenna gain and the total 
power radiated. EIRP is determined as [Ref. 4: p. 62]: 

EIRP = P rad D 0 (W) (4.11) 

Assuming that the power delivered to planar arrays 
constructed of M x N isotropic radiating elements is 
equivalent to the radiated power, we estimate the radiation 
resistance (R r ) of the array as [Ref. 3: p. 55]: 

s, - “f <fi> (4.12) 

Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a planar array (R in ) is 
assumed to be equivalent to the radiation resistance: 

R in = R z (fi) (4.13) 

The gain (G) of a planar array cannot be calculated given 
the information assumed to be available, but an idealized 
value can be derived. Another factor affecting the gain is 
the assumption made when determining the total antenna 
efficiency (e t ) . The total antenna efficiency for an antenna 



23 



is [Ref. 3: pp. 44-45]: 



c c = c r e cd = l - | r | 2 ( dimensionless ) (4.14) 



where 



e cd = 1 (dimensionless) 



(4.15) 



e r = 1 - I r I 2 ( dimensionless ) (4.16) 

r = — — (dimensionless) (4.17) 

R in + Z o 

In (4.14), e cd is the conduction and dielectric efficiency, 

which is unity for an ideal lossless antenna. In (4.17), Z G 

is the characteristic impedance of the transmission line and 
r is the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(4.15) when ideal lossless antenna is assumed. The ideal gain 
of a planar array in dimensionless and decibel quantities is 
then [Ref. 3: pp. 43-44]: 

G = z t D 0 ( dimensionless ) (4.18) 



G(dB) = 10log 10 ( e c D 0 ) (dB) (4.19) 

The maximum effective aperture (A ero ) for a planar array 
is derived by assuming polarization mismatches are negligible 
(this application is constructed using isotropic radiating 
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elements for which polarization loss factor (PLF) and antenna 
polarization are not calculated) [Ref. 3: p. 63]: 




( 4 . 20 ) 



PLF = 1 ( dimensionless ) 



(4.21) 



The effective height (h em ) of the planar array is now 
estimated by [Ref. 6: p. 42]: 



In (4.22), r) 0 is the characteristic impedance of free space 

(120tt) . 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is in the x-z plane; 
the second is in the y-z plane; the third display generated is 
a plane perpendicular to the x-y plane when 0 = 7r/4 ; the last 
display is a plane perpendicular to the x-y plane when 0 = <p a . 
The last display is redundant when 0 is either 0, tt/ 4, or n/2. 

Results obtained with this application are the same as 
those found in [Ref. 3: pp. 260-274]. Directivity was 

calculated using the approximate equation (4.10). An 
alternative approximation for directivity is given by either 
[Ref. 3: p. 272-273]: 




(4.22) 
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D a = tc COS ( 0 O ) D x D y 



( dimensi onl ess) 



(4.23) 



or 

2 

D “ ( dimensionless ) (4.24) 

° Q A (rads 2 ) 

In (4.23), D x and D y are the directivities of broadside linear 
arrays of length and number of elements Md x , M, and Nd y , N in 
the x and y directions, respectively. In (4.24), fi A is the 
beam solid angle. 

Table 4.1 is a comparison of directivity computed for a 
10x10 planar array. The reference directivity is derived 
using (4.24) while the calculated directivity is from (4.10). 



Table 4.1 Planar Array Data Comparison 




Figure 4.2 and Figure 4.3 are comparisons of array factor 
patterns for a 5x5 planar array with patterns from [Ref. 3] 
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and the calculated patterns. Figure 4.2 is a representation 
of the array factor patterns in the x-z and y-z planes, and 
the scale of this figure has been increased in order to 
increase the clarity of the pattern. Figure 4.3 is the array 
pattern in the plane at 0 = 45°, where the maximum beam is 
directed at <p 0 = 45°. The array patterns produced by this 
application are identical to the patterns given in [Ref. 3]. 
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FIGURE 4.2 Array factor pattern of a 
excited, planar array with d x = d y = 
45°, in the x-z and y-z planes 



5x5 element, equally 
A/2 m, 6 0 = 30° ,<p 0 = 
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FIGURE 4.3 Array factor pattern of a 
excited, planar array with d x = d y = 
45' , in the plane 0 = 45' 



5x5 element, equally 
X/2 m, 0 O = 30', 0 O = 
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V. THE CIRCULAR ARRAY ANTENNA 



Circular arrays, or ring arrays, are antennas in which 
the radiating elements are placed on a circle with no elements 
positioned inside the circle. Circular antenna arrays with 
elements positioned inside the circle are planar arrays, and 
these antennas are addressed in Chapter IV. Circular array 
antenna are used for radio direction finding, air and space 
navigation, underground propagation, radar, and sonar systems 
[Ref. 3: p. 274]. As with linear and planar arrays, circular 
arrays with uniformly spaced, equally excited elements are 
narrowband antenna systems due to the dependance of the number 
of elements (N) and the circumference of the circle to the 
wavelength (1) [Ref. 3: p. 277]. 

The following inputs are required to analyze a circular 
array : 

N = number of isotropic radiating elements 
a = radius of the circle 
f = frequency of interest * 

0 O = direction of main lobe at 0 = 0 O * 

0 O = direction of main lobe at cp = <p Q * 

I 0 = antenna feed current 

Z 0 = characteristic feed impedance " 

The first two inputs are data primarily obtained from visual 
data such as photographs, and the * indicates input parameters 
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that are either known or estimated. The planar array geometry 
is illustrated in Figure 5.1. 




FIGURE 5.1 Geometry of N-element circular array of 
isotropic point sources 
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This Mathcad application is written for N-element 
circular arrays with uniformly spaced, equally excited, 



isotropic radiating elements on the x-y plane. The 



polarization of the array is not calculated since isotropic 
point sources are used, and the antenna input resistance (R in ) 
is assumed to be equivalent to the radiation resistance (R r ) . 
This assumption results in the computation of ideal gain (G) . 
As with the linear and planar arrays, the radiation pattern of 
a specific element type is neglected in the analysis of the 
circular array since in normal usage it will have little 
effect for a circular array with a large number of elements. 
However, as with both the linear and planar arrays, the user 
of this application should be able to use information 
regarding the type of element and the orientation of the 
elements to narrow the range of possible 0 o 's and <p Q f s. 

The normalized array factor (AF) is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [Ref. 3: pp. 
274-278; Ref. 4: 350-354]: 




j [ka sin6 cos aj 



( dimensionless ) (5.1) 



where 




(5.2) 



a n - sin0 o cos ( 4> 0 - <I> n ) ( radians ) 



(5.3) 
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Equation (5.1) is the normalized array factor. Equation (5.2) 
is the angular position of the nth element on the x-y plane. 
Equation (5.3) is the phase excitation of the nth element with 



the position of the main beam directed at 0 = 0 O and <p = <p ol 
and k is the wavenumber ( 2n/X ) and a is the radius of the 
array . 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 



The minimum distance to the far-field is found by taking the 
maximum of (5.4), (5.4), and (5.6), where the diameter of the 
circle is taken to be the maximum dimension of the antenna. 

To determine the directivity (D 0 ) for a circular array, 
the approximate equations of radiation intensity (U(0,0)) for 
arrays and radiated power (P rad ) are used. Radiation intensity 
is [Ref. 3: pp. 229-233]: 



r ^ 1.6 X ( meters ) 



(5.4) 



r £ 5 ( 2 a ) ( meters ) 



(5.5) 



i 2: 2J_2a) ( meters ) 

A 



(5.6) 



t/(0,4>) = ( AF ) 2 ( W/ solid angle) 



(5.7) 



The radiated power is [Ref. 3: pp . 28]: 




(5.8) 
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In (5.8), P rad is computed by integrating over all angles above 
the x-y plane (0 < 0 < n/2) [Ref. 7: p. 78]. This assumes 

that the array has a pattern maxima above the x-y plane with 
no radiation below the x-y plane. The directivity for the 
circular array is [Ref. 3: pp. 229-233]: 

D = - - ^ / ' 7iax ( dimensionless ) (5.9) 

P iad 

In (5.9), U raax is unity and occurs at 0 = 0 O and <p = <p D since 
the array factor is normalized. 

Effective isotropic radiated power (EIRP) is determined 
by [Ref. 4 : p. 62 ] : 

EIRP = P zad D 0 (W) (5.10) 

Assuming that the power delivered to a circular array of 
N isotropic, uniformly spaced radiating elements is equivalent 
to the radiated power, we can estimate then the radiation 
resistance (R r ) of the array as [Ref. 3: p. 55]: 

R r = 2 I (ft) (5.11) 

l-U 

Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a circular array (R in ) is 
estimated to be equivalent to the radiation resistance: 

R in = R z (ft) (5.12) 
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Since the input resistance is estimated, the gain (G) of 
a circular array cannot be calculated, but an ideal value can 
be derived. As with linear and planar arrays, addressed in 
Chapter III and Chapter IV, respectively, the total antenna 
efficiency (e t ) is [Ref. 3: pp. 44-45]: 

e t = e r e cd = l-|r| 2 ( dimensionless ) (5.13) 



where 



e cd = 1 ( dimensionless ) 



(5.14) 



e r = l-|r| 2 (dimensionless) (5.15) 

r = — ^ — £2 ( dimensionless ) (5.16) 

Rin + Zo 

In (5.14), e cd is the conduction and dielectric efficiency, 
which is unity for an ideal lossless antenna. In (5.16), Z 0 
is the characteristic impedance of a transmission line, and r 
is the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(5.15) when an ideal lossless antenna is assumed, and the gain 
of a circular array is [Ref. 3: pp. 43-44]: 

G = t c D 0 ( dimensionless ) (5.17) 



G(dB) = 10log 10 ( z t D a ) ( dB ) (5.18) 
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The maximum effective aperture (A em ) of a circular array 
is derived assuming polarization mismatches are negligible 
(this Mathcad application is constructed using isotropic 
radiating elements for which polarization loss factor (PLF) 
and antenna polarization are not calculated) [Ref. 3: p. 63]: 

A = 91121 PLF (m 2 ) (5.19) 

em 4 71 

PLF = 1 ( dimensionless ) (5.20) 



The effective height (h e J of a circular array is 
estimated by [Ref. 6: p. 42]: 






*r A em 



(m) 



(5.21) 



In (5.21), ho is the characteristic impedance of free space 
(12077) . 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is the x-z plane; the 
second is the y-z plane; the third is a plane perpendicular to 
the x-y plane when <p = 77/4; the last is a plane perpendicular 
to the x-y plane when 0 = 0 O . The last display is redundant 
when 0 is either 0, 77/4, or 77/2. 

The directivity of a uniform circular array approaches 
the number of elements (N) as the radius (a) of the array 
becomes very large [Ref. 3: p. 277]. Table 5.1 provides the 
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directivity of a uniform circular array as a function of the 
radius of the array. 



Table 5.1 Circular Array Directivity Data 



RADIUS (m) 


NUMBER OF 


DIRECTIVITY (D 0 ) 




ELEMENTS 




1.59 


10 


13.70 dB 


10 


10 


12.83 dB 


20 


10 


12.89 dB 


50 


10 


11.14 dB 


100 


10 


10.33 dB 



Figures 5.2 and 5.3 provide comparisons of the array 
radiation patterns computed for a 10-element, uniformly 
spaced, equally excited circular array with a radius of 1.59 
meters with those in [Ref. 3]. Figure 5.2 is the radiation 
pattern in the x-z plane, and the radiation pattern in the y-z 
plane is shown in Figure 5.3. For these patterns, 0 O = <p D = 
0°. The array radiation patterns generated by this 
application are identical to the patterns given in [Ref. 3]. 
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-1 - 0,6 - 0.6 - 0,4 - 0,2 0 0,2 0.4 0,6 0,8 1 



FIGURE 5.2 Array factor pattern of a 10-element, uniformly 
spaced, equally excited circular array with a = 1.59 m, 
and 0 O = 0 O = 0°, in the x-z plane 
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FIGURE 5.3 Array factor pattern of a 10-element, uniform 
spaced, equally excited circular array with a = 1.59 m, 
and 0 O = <p Q - 0°, in the y-z plane 
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VI. THE FOLDED DIPOLE ANTENNA 



The folded dipole is a center-fed wire antenna that 
performs in a manner similar to that of a linear, center-fed 
thin-wire antenna with the added advantage of a more easily 
matched input impedance [Ref. 8: p. 7-36]. The folded dipole 
antenna consists of two parallel dipoles connected at the ends 
to form a narrow wire loop. The folded dipole antenna is 
illustrated in* Figure 6.1. 




FIGURE 6.1 Geometry of folded dipole 
antenna 
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The folded dipole can be assumed to behave as a thin-wire 
dipole when the following are satisfied: the length of the 
antenna is much greater than the diameter of the feed 
conductor (L > 2a); the distance between the conductors is 
greater than the radius of the feed conductor (d > a) ; 
(27rd/A.) 2 < 1; all conductors are assumed to be perfect [Ref. 
9: p. 172]; 

This Mathcad application calculates the parameters of a 
folded dipole by assuming the antenna behaves as a thin-wire 
dipole with an equivalent radius (a e ) . The current maximum is 
normalized to one amp. The equivalent radius is [Ref. 8: p. 
7-38] : 



a 



e 



exp 



a 2 In (a) + b 2 ln(b) + 2abln(d) 

(a + b ) 2 



( meters ) (6.1) 



This application provides the characteristics of the antenna 
both in free space and positioned horizontally over the earth. 
Although the folded dipole positioned vertically over the 
earth occurs frequently in VHF/UHF communications, it is not 
addressed in this research. 

The following inputs are required to analysis a folded 
dipole : 

a = radius of feed conductor 
b = radius of second conductor 
d = spacing between conductors 
L = length of antenna 
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h = height of antenna above ground plane 
f = frequency of interest 

e cd = conduction/dielectric efficiency of conductors 
e r = relative dielectric constant of ground plane 
a = conductivity of ground plane 

r ff = distance for far-field parameter calculations 
a w = incoming wave electric field unit vector for antenna 
in free space 

cr a = unit polarization vector for antenna in free space" 
crh w = incoming wave electric field unit vector for 
horizontally positioned antenna 
0h r = antenna polarization direction for horizontally 
positioned antenna 

The first five inputs are physical dimensions obtained through 
photographs, the " indicates input parameters that are either 
known or estimated. 

A. FOLDED DIPOLE IN FREE SPACE 

Assuming the folded dipole behaves as a thin-wire dipole 
with a e < d < L, we have the far-field electric field 
intensity of the antenna aligned with the z-axis as [Ref. 3: 



In (6.2), r) 0 is the intrinsic impedance of free space, equal 



p. 120]: 
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to 1207T, and k = 2 n/X is the wavenumber. The normalized 

electric field is found by dividing (6.2) by the maximum value 
of the electric field which occurs at 6 = n/2 for L less than 
about 1.25A. Clearly, the radiation pattern of the folded 
dipole is identical to that of a conventional dipole of length 
L. 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 

r ^ 1 .6 A ( meters ) (6«3) 

r £ 5 L {meters) (6-4) 



r ^ 



2 L 2 
1 



( meters ) 



(6.5) 



The minimum distance to the far-field is found by comparing 
the values of (6.3), (6.4), and (6.5) and selecting the 
maximum value. The distance between the conductors is not 
taken into account in determining the minimum distance to the 
far-field since it is assumed that the folded dipole has an 
over all effective radius (a e ) . 

Directivity (D 0 ) for the folded dipole is found by 
determining the radiation intensity (U(0)) and radiated power 
(P rad ) . Since the radiation pattern is symmetric in the x-y 
plane, U(0) and P rad are functions of 0 only. Radiation 
intensity is [Ref. 3: p. 28]: 
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U(Q) 



( 6 . 6 ) 




The radiated power is [Ref. 3: pp. 28]: 



P rad = 2nf n U(d) sin (6) dd (W) 
J 0 



(6.7) 



Directivity for a folded dipole is [Ref. 3: pp. 29]: 



4 Tt U. 



max 



( dimensi onl ess) 



( 6 . 8 ) 



D, 



O 



where U max = U(tt/2) . 

Effective isotropic radiated power (EIRP) is determined 
by [Ref . 4 : p . 62 ] : 



The folded dipole is unique in its capability to act as 
a step up impedance transformer when the antenna is resonant. 
The antenna effectively operates as an unbalanced transmission 
line, and it can be analyzed by decomposing the current into 
two distinct modes: a transmission line mode [Figure 6.2] and 
an antenna mode [Figure 6.3]. [Ref. 3: p. 341; Ref. 4: pp. 
205-207] 



EIRP = P rad D 0 (W) 



(6.9) 
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V 

2 



FIGURE 6.2 Transmission 
line mode 




The input impedance (Z in ) at the feed for the folded 
dipole is given in (6.10), and the radiation resistance (R in ) 
for the antenna is the real component of (Z in )) [Ref. 3: 
p.342] : 

4Z.Z, 



Z, 



‘t^d 



in 2 Z d + Z c 



(q; 



( 6 . 10 ) 
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In (6.10), Z t is the input impedance of the transmission line 
mode and is obtained from the impedance transfer equation 



Z c = J2 0 tan|-^) (Q) 



( 6 . 11 ) 



where Z Q in (6.11) is the characteristic impedance of a two- 
wire transmission line and is given in (6.12) [Ref. 3: pp. 

341-342; Ref. 8: p. 7-39]: 

Z c = — cosh' 1 / — ) (Q) (6.12) 

71 \2/aB) 

The input impedance of a linear dipole of length L and 
equivalent radius a e is calculated by [Ref. 3: pp. 124, 127, 

294 , and 342 ] : 

Z d = R d + jX d (Q) (6.13) 



R d = 



sm- 



l k_L \ 



(Q) 



(6.14) 



= 



X r 



sm* 



kL\ 

2 



(Cl) 



(6.15) 



R 



r = y|[Y + ln(icL) - CjUL)] 



2 TI 



sin{ kL) [S i (2kL) -2S i {kL)} 



(6.16) 



2;t 



cos ( kL) |y + + C i ( 2 kL) -2 C i ( kL) 



46 



(6.17) 



X r = — [2 S, ( kL) +cos ( kL) [2 5, (kL) - S i ( 2 k L) ] 
1 4 7t l 1 1 



4 n 



sin ( kL) 



2C 2 (kL) -C 1 (2kL) -C 1 



2 ka e 
L~ 



In (6.16) and (6.17), C x (x) and S^x) are the cosine and sine 
integrals, respectively, [Ref. 3: pp. 124 and 743-746]: 



C 2 (x) 



[ x cos(jl dT = Y+ln(x) 

t hi 



( -1) ft-J-x 2 ” 
2 n ( 2 n ) ! 



(6.18) 



w - ■ t 



( -D^X 12 "' 11 
(2 n - 1 ) (271-1) ! 



(6.19) 



and y = 0.57721 (Euler's constant). 

Mathcad computes summations more quickly than integrals. 
In this application, the summation forms of (6.18) and (6.19) 
are used with the index n ranging from 1 to 50 to ensure 
proper convergence. The maximum value of the argument of 
either C^x) or Sj(x) is x = 4 ttL/A., and for L = X the 
succeeding terms in the summation rapidly decrease for n > 20. 

The gain (G) of a folded dipole is obtained as the 
product of the antenna reflection efficiency (e L ) and the 
directivity of (6.8) where the antenna reflection efficiency 
is [Ref. 3: pp. 43-45]: 

C c = e r e cd = e i - | F | 2 ) ( dimensionless ) (6.20) 
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where 



c z = i - | r | 2 ( dimensionless ) (6.21) 

p = 3m — £2 ( dimensionless ) (6.22) 

Rin + Zo 

In (6.21) e r is the mismatch efficiency and e cd is the 
conduction and dielectric efficiency in (6.22). In (6.22), Z Q 
is the characteristic impedance of the transmission line given 
in (6.12). Therefore, the gain of the antenna is: 

G = t c D 0 ( dimensionless ) (6.23) 

G(dB) = lOlog 10 (t c D o ) ( dB ) (6.24) 

The folded dipole in free space aligned on the z-axis is 
horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave's unit 
polarization vector (a w ) , the polarization loss factor (PLF) 
is [Ref. 3: pp. 48-53]: 

PLF = \a w - o* a \ 2 ( dimensionless ) (6.25) 

The maximum effective aperture (A em ) of a folded dipole 
is [Ref. 3 : p. 63] : 

A em = G [ X) 2 PLF (in 2 ) (6.26) 

4 K 
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The effective height (h em ) of the folded dipole is [Ref. 



6 : p. 42 ] : 



= 2 



Rz A en, 



\ 



(m) 



(6.27) 



The bandwidth of the folded dipole antenna is dependant 
on the length of the wires and has better bandwidth 
characteristics than a conventional dipole of the same length. 
It can be assumed that the bandwidth of the folded dipole is 
essentially the same as that of a conventional dipole of an 
radius a e (a < a e < d/2) [Ref. 3: p. 346]. The bandwidth is 
increased when the distance between the conductors (d) is 
increased and/or when the conductor radii a and b are 
increased [Ref. 10: p. 190]. 

Tables 6.1 and 6.2 are comparisons of measured and 
calculated input impedance data for the folded dipole as a 
function of its electrical length. The radii of the 
conductors (a, b) are the same (0.0005 m) while the spacing 
between the conductors (d) is held constant at 0.00625 m and 
X = 1 m [Ref. 8: pp. 7-40 - 7-41]. 
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TABLE 6.1 Folded Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


CALCULATED 


LENGTH (L) 


Re ( z in ) n 


Re ( z in ) n 


X/n 


2000 


6533 ** 


2X/n 


2400 


2592 


3 X/n 


H • 0 


0.1 


4X/ir 


300 


324 


5X/TT 


1700 


1685 


6X/TT 


£•• 


1»0 


1 X/n 


100 


100 
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TABLE 6.2 Folded Dipole Antenna Data Comparison 



ANTENNA 
LENGTH (L) 


MEASURED 
Im(Z in ) n 


CALCULATED 
Im(Z in ) n 


X/ir 


5000 


5437 


2 A/ 7T 


-650 


-1204 ** 


2X/n 


-90 


-85 


AX/n 


1000 


1061 


5X/ir 


700 


793 


&X/1T 


-190 


-174 


iX/n 


600 


633 



To further test the validity of the Mathcad application, 
a numerical antenna analysis program was also used. The 
program used was ELNEC, a powerful, easy-to-use program for 
modeling and analyzing a wide variety of antenna types 
including ground effects and parasitic structures. The 
fundamental computation portion of ELNEC is the same as 
MININEC (Version 3), that was developed by the Naval Ocean 
Systems Center [Ref. 11: pp. 5-7]. 

Table 6.3 is a comparison of the gain for a folded dipole 
in free space with different dimensions [Ref. 8: p. 3-29]. 
For each case, the radii of the conductors (0.0005 m) and the 
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distance between the conductors (12.5a m) are the same. As 
can be seen, the comparison between the gains obtained from 
various sources is very good. 



TABLE 6.3 Folded Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


ELNEC 


Mathcad 


DIMENSION 


Gain (G) 


Gain (G) 


Gain (G) 


L = A/4 m 


1.71 dB 


1.70 dB 


1.64 dB 


d = L/13 m 








L = A/2 m 


2.14 dB 


2.20 dB 


2.15 dB 


d = L/ 2 5 . 5 m 


j 







The computed parameters in Tables 6.1, 6.2, and 6.3 are 
essentially the same as those measured or computed with ELNEC 
with exception of the ones marked with **. These discrepancies 
result in using the transmission line model and equivalent 
radius to analysis the folded dipole in this application. 

As previously mentioned, the folded dipole is a center- 
fed wire antenna that performs in a manner similar to a 
linear, center-fed, thin-wire antenna. Therefore, the 
radiation pattern of a folded dipole is also similar to that 
of a thin-wire antenna. The folded dipole radiation pattern 
in free space is illustrated in Figure 6.4. 
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FIGURE 6.4 Radiation Pattern folded dipole in free space 
with L = A/2 m, a = b = 0.0005 m, and a e = 0.004 m 



Figure 6.5 is a comparison of the radiation pattern of a 
folded dipole in free space obtained using ELNEC and the 
Mathcad application. The patterns are not identical, but are 
similar and oriented in the same direction as expected. The 
difference can be attributed to the fact that ELNEC takes into 
account mutual coupling between the conductors while the 
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Mathcad application models the folded dipole 
center-fed dipole. 



as 



a 



thin-wire 




Figure 6.5 Comparison of radiation patterns using ELNEC 
and Mathcad of a Folded dipole in free space with L = A/2 
m, a = b = 0.0005 m, and a e = 0.004 m 
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B. FOLDED DIPOLE POSITIONED HORIZONTALLY OVER THE EARTH 

Through the use of image theory, the far-field electric 
field intensity of a folded dipole oriented horizontally above 
a flat earth can be obtained. The 6 component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 
space is [Ref. 3: pp. 143-144]: 



From image theory, the 0 component of the electric field 
over ground is [Ref. 4: pp. 229-235]: 



In (6.28) and (6.29), E e (0,<£) is valid only above the ground 
(0 < 0 < n/2) [Ref. 3: pp. 135-142]. In (6.29), h is the 
distance from the earth to the antenna and the vertical 
reflection coefficient (T v ) of the ground is [Ref. 4: pp. 231- 
233 ] : 






v^l - sin 2 (0) sin 2 (4>) 



£0 gp (0) =Ed (6,4>)(e jkhcos(e ‘ - r v e- jkhcos,e) ) (V/m) (6.29) 



e / r cos(0) - s ft r I - sin 2 ( 0 ) 



(dimensionless) (6.30) 



e'j.cos ( 0 ) + r - sin 2 ( 0 ) 



where 



I 




O 



( dimensionless) 



(6.31) 
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In (6.30) and (6.31), (e' r ) is the relative complex effective 
dielectric constant of the ground, (e r ) is the relative 
dielectric constant of the ground, and (a) is the conductivity 
of the ground. For a perfect ground plane, a -* °o and r v = 1. 

The # component of the electric field of a horizontal 
dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 
135-142] : 



I e~ jkz “ 

£•<!> ( 0 , 4 >) =jr\ ° 



2 nr 



ft 



cosf^sin (0) sin(4>)J -cosj^j 



\fl - sin 2 (0) sin 2 (4>) 



(v/m) (6 - 32) 



Now from image theory, the 4> component of the electric field 
over earth is [Ref. 4: pp. 229-235]: 

£<*>^(0) =£’<M0,4>)(e^ cos(!,) + r\e-^ ICO£<0) ) (v/m) (6.33) 



In (6.33), (T h ) is the horizontal reflection coefficient of 
the ground and is [Ref. 4: pp. 231-232]: 



cos ( 0 ) - Je,' - sin 2 ( 0 ) 

( dimensionless ) (6.34) 

cos(0) + ^e 7 r - sin 2 ( 0 ) 



When o -» oo, r h = -1. 

The radiation intensity (U h (0,0)) of a folded dipole over 
a earth must take into account both the 0 and components of 
the electric field. Therefore, the radiation intensity is 
[Ref. 3: p. 28]: 
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U h (6) = ^-\\E8 gp {Q) | 2 + |£’4> 5p (e) | 2 ] ( W/ solid angle) (6.35) 

The radiated power (Ph rad ) follows (6.7) except that the field 
is present only above the earth. Thus, 

71 

Ph rad = 2 tt f 2 U h ( 0 ) s in ( 0 ) dQ (W) (6.36) 

J o 

Directivity (Dh 0 ) and effective isotropic radiated power 
(EIRPh) for a horizontal folded dipole over the earth are 
identical to (6.8) and (6.9), respectively, with the results 
of (6.36) inserted. 

Input resistance (Rh in ) for the folded dipole over the 
earth plane is found by (where I c in (6.32) is normalized to 
one amp) [Ref. 3: p. 124]: 




where the radiation resistance (Rh r ) is, 

2 Ph 

Rh r = -—^ (Q) (6.38) 

\lo\ 

The voltage reflection coefficient (TH) , gain (Gh) , 
maximum effective aperture (Ah em ) , and maximum effective height 
(hh em ) for a horizontal folded dipole over the earth are 
determined using the equations for those parameters for a 
folded dipole in free space, as listed above. 
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The Polarization Loss Factor ( PLFh) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (0h p ) is estimated. The antenna 
polarization vector (oh a ) is then calculated, and estimating 
the unit polarization vector of the incoming wave (ahj , the 
PLF is found from (6.25). 

Table 6.4 is a comparison of measured [Ref. 8: p. 3-29], 
ELNEC, and Mathcad software gains for a horizontal folded 
dipole over a reflecting sheet. The radii of the conductors 
(0.0067 m) are the same. As can be seen, the Mathcad software 
gain is similar to the measured and ELNEC gains. 



TABLE 6.4 Folded Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


ELNEC 


Mathcad 


DIMENSION 


Gain (G) 


Gain (G) 


Gain ( G ) 


L - A/2 m 








d = L/25.5 m 


7.14 dB 


6.98 dB 


7.19 dB 


h = A/8 m 


L 







No input impedance data has been found to compare a 
folded dipole oriented horizontally over the earth with that 
computed by this Mathcad application. 

Figures 6.6 and 6.7 are the H-plane and E-plane radiation 
patterns of a folded dipole oriented horizontally over the 
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earth, respectively. Since the folded dipole is similar to a 
thin-wire center-fed dipole of length L, the illustrations 
below are similar to those of a dipole of length L. 




FIGURE 6.6 H-plane radiation pattern of a folded dipole 
with L = X/2 m, a = b = 0.0005 m, d = 0.00625 m, h = 0.5 
m, c = 0.01 S/m, e r = 15 
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FIGURE 6.7 E-plane radiation pattern of a folded dipole 
with L = A/2 m, a = b = 0.0005 m, d = 0.00625 m, h = 0.5 
m, o = 0.01 S/m, e r = 15 



Figures 6.8 and 6.9 are comparisons of the H-plane and E- 
plane radiation patterns between those obtained with ELNEC and 
Mathcad software of a folded dipole oriented horizontally over 
the earth, respectively. As expected, the patterns are 
similar. As previosly mentioned, the difference can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the conductors while the Mathcad application 
is modeled after a thin-wire, center-fed dipole. 
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Figure 6.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a folded dipole with L = A/2 m, a = 
b = 0.0005 m, d = 0.00625 m, h = 0.5 ra, a = 0.01 S/m, e r = 
15 




Figure 6.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a folded dipole with L = A/2 m, a = 
b = 0.0005 m, d = 0.00625 m, h = 0.5 m, a = 0.01 S/m, e r = 
15 
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VII. THE CAGED DIPOLE ANTENNA 



The caged dipole antenna is used in communications for 
its broadband antenna characteristics. Its performance is 
similar to that of a thick cylindrical dipole. Since thick 
cylindrical dipoles are usually too clumsy to use except at 
very short wavelengths, thick dipoles are replaced with a 
caged dipole which consists of a number of thin wire dipoles 
arranged to form a circular cross section. The antenna is 
usually supported by a central conductor and support members 
that provide rigidity but contribute little to the electric 
properties of the caged dipole. The caged dipole (without the 
central conductor and the support members) is illustrated in 
Figure 7.1 [Ref. 12: pp. 172-174]. 
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This Mathcad application models the caged dipole antenna 
as a circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles. This approach is taken to account 
for the number of conductors that are present in the caged 
dipole. As with the arrays previously considered, mutual 
coupling between the dipoles is neglected. In this manner, 
this application simulates the characteristics of a thick 
cylindrical dipole. The thin-wire dipoles are modeled after 
finite length dipoles in free space that have their length 
much greater than their diameter (L > wire diameter) , and the 
dipoles are assumed to be perfect conductors. Since the 
central conductor and the support wires have little effect on 
the electrical properties of a caged dipole, mutual coupling 
from these components are not taken into account in this 
application. As seen in Figure 7.1, the inner ends of a caged 
dipole, the feed points, maybe coned in order to reduce base 
capacity [Ref. 12: pp. 172-174]. This effect is also not 
addressed . 

This application provides the characteristics of the 
antenna in free space as well as and oriented both vertically 
and horizontally over the earth. 



A. CAGED DIPOLE IN FREE SPACE 

The following inputs are required to analyze a caged 
dipole in free space: 
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N = number of conductors 



a = radius of antenna 
L = length of antenna 
f = frequency of interest 
I c = antenna feed current 
Z 0 = input impedance 

e cd = conduction/dielectric efficiency of conductors 
r ££ = distance for far-field parameter calculations 
a w = incoming wave electric field unit vector for 
antenna 

a a = unit polarization vector for antenna 
The first three inputs are physical dimensions obtained 
primarily from photographs, and the * indicates input 
parameters that are either known or estimated. 

Assuming the caged dipole behaves as a circular array 
with equally excited (I n ), uniformly spaced elements in the x- 
y plane with a radius of (a) , we have the array factor 
( AF (6,0) ) [Ref. 3: pp . 274-278]: 

AF(0,4>) = J n e j I * asin0cos<<l ’~ ,I> " )] ( dimensionless ) (7*1) 

n- 1 

where the angular position in the nth element in the x-y plane 
is : 

n = 1,2, . . . N (radians) (7.2) 
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Equation (7.1) is derived from: 



AF(0,<t>) = ^ I n e j [Jcasin6cos(<t ’~' t> « ) *°" ] ( dimensionless ) (7.3) 

n* 1 

where the phase excitation of the nth element is: 

<x n = 0 = -kasin0 o cos (4> 0 -«> n ) ( radians ) (7.4) 



In (7.3), and the phase excitation of the nth element (a n ) is 
zero as the result of desiring the main beam of the array to 
be directed perpendicular to the length of the antenna (0 O = 
7r/2). In (7.1), k is the wavenumber (2ir/X). 

Assuming each conductor of a caged dipole behaves as a 
thin-wire dipole, we have the far-field electric field 
intensity of a single conductor aligned with the z-axis as 
[Ref. 3 : p. 120 ] : 



F(0) = JTj J " e 



-jkif 



2 nr 



ft 



cos| cos (0) j - c osj-y 
sin (0) 



(V/m) 



(7.5) 



In (7.5), r ] o = 1 2 Off and is the intrinsic impedance of free 
space . 

As previously discussed in Chapter III, the total 
electric field of an array constructed with identical elements 
is the product of the field of a single element and the array 
factor as a consequence of the principal of pattern 
multiplication [Ref. 3: p. 207]. Therefore, the total 



65 



electric field of a caged dipole in free space is approximated 
as : 

£ c <e,4>) = £(6) AF(0,<t>) (V/m) (7.6) 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 

r s 1.6 A. ( meters ) (7.7) 



r ;> 5\/L 2 + (2a) 2 



r £ 



2 (L 2 + (2a) 2 ) 

X 



(meters) 



(meters) 



(7.8) 



(7.9) 



The minimum distance to the far-field is found by taking the 
maximum of (7.7) through (7.9). In (7.8) and (7.9), the 
maximum dimension of the antenna is the diagonal length of the 
array . 

Directivity (D 0 ) for a caged dipole is found by 
determining the radiation intensity (U(0,0)) and radiated 
power (P rad ) . Radiation intensity is [Ref. 3: p. 28]: 

t/(0,4>) = Ift |£ c ( 0 , <{> ) | 2 (w/ solid angle) (7.10) 

2 'Ho 

The radiated power is [Ref. 3: pp. 28]: 

p rad = f 2n f n U(e,4>) sin(0) c/ed(p ( IV) (7.11) 

Jo Jo 
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Directivity for a caged dipole is [Ref. 3: pp. 29]: 



D a = 4n U ™ ax ( dimensionless ) (7.12) 

^rad 

where U ffiax = U(tt/2 , 0) . 

Effective isotropic radiated power (EIRP) is determined 
by [Ref. 4 : p. 62 ] : 

EIRP - P iad D 0 (IV) (7.13) 



Input resistance (R in ) for a caged dipole in free space 
is found from [Ref. 3: p. 124]: 

R, 



R, 



sirr 



2 , f) 



(Q) 



(7.14) 



where the radiation resistance (R r ) is: 

2 P r 



*r = 



rad 



UJ 2 



(Q) 



(7.15) 



The gain (G) of a caged dipole is obtained as the product 
of the antenna reflection efficiency (c L ) and the directivity, 
where the antenna reflection efficiency is ([Ref. 3: pp. 43- 
45] : 

e t = = e«,(l-|IT) ( dimensionless) (7.16) 



where 

e r = l-|r | 2 ( dimensionless ) (7.17) 
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( dimensionless) 



(7.18) 




The parameter e r is the mismatch efficiency, e cd is the 
conduction and dielectric efficiency which is unity for a 
lossless conductor, and Z c is the characteristic impedance of 
a transmission line. Therefore, the gain of the antenna is: 



A caged dipole in free space aligned parallel to the z- 
axis is horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave's unit 
polarization vector (a w ) , the polarization loss factor (PLF) 
is [Ref. 3: pp. 48-53]: 



The maximum effective aperture (A em ) of a caged dipole is 
[Ref. 3: p.63]: 



The effective height (h em ) of the caged dipole is [Ref. 
6: p. 42]: 



G = z t D a ( dimensionless ) 



(7.19) 



G(dB) = 10 log 10 (e c D a ) ( dB ) 



(7.20) 



PLF = | 5 




( dimensionless ) 



(7.21) 




(7.22) 



h 




(7.23) 
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The bandwidth of a thick cylindrical dipole antenna is 
dependent on the ratio of the length of the antenna to the 
diameter of the antenna (L/2a) . When L/2a ~ 5000 the 
acceptable bandwidth of a thick cylindrical dipole is about 3% 
of the center frequency, and when L/2a ~ 260 the bandwidth is 
about 30%. Hence, if we assume a linear relationship between 
bandwidth and L/2a, and since the caged dipole's performance 
is similar to that of a thick cylindrical dipole, then the 
bandwidth of a caged dipole is [Ref. 3: p. 333]: 

BW = W - flow (Hz) (7.24) 

where the upper and lower frequencies are: 

W = fc + F (j^) f c < 7 - 25 ) 

fio» = fc - (Hz) ( 7 * 26 ) 

In (7.25) and (7.25), f c is the center frequency and F(L/2a) 
is the linear relationship function between antenna length and 
diameter. This relationship is obtained by solving y = mx + 
f c for m using Xj = 5000, y x = 3% and x 2 = 260, y 2 = 30%. The 
relationship becomes: 

F (la) = _5 ‘ 696 ' 10 ' 5 (' 2a) + °' 3148 ( dimensionless) (7.27) 
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Table 7.1 is a comparison of the gain of a thick 
cylindrical dipole with that of a caged dipole (N = 8) in free 
space computed both by ELNEC and Mathcad software [Ref. 8: p. 
2-38]. As previously mentioned in Chapter VII, ELNEC is an 
easy-to-use numerical analysis program for modeling and 
analyzing different antenna types [Ref. 11: pp. 5-7]. As can 
be seen, the comparison between the measured, ELNEC, and 
Mathcad application values are similar. 



TABLE 7.1 Caged Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


ELNEC 


Mathcad 


DIMENSIONS 


Gain (G) 


Gain (G) 


Gain (G) 


L = k/2 m 








a = 0.025 m 


2.14 dB 


1.87 dB 


1.88 dB 


N = 8 









Figure 7.2 is the radiation pattern of a caged dipole (N 
= 8) in free space. It is essentially equivalent to that of 
a center-fed, thick cylindrical dipole having the same 
conductor length and diameter [Ref. 8: p. 2-38]: 
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FIGURE 7.2 Radiation pattern of a caged dipole in free 
space with L = A/2 m, N = 8, and a = 0.125 m 



Figure 7.3 is a comparison of the radiation pattern of a 
caged dipole (N = 8) obtained using ELNEC and Mathcad. The 
patterns are not identical, but are similar and oriented in 
the same direction as expected. The difference can be 
attributed to the fact that ELNEC takes into account mutual 
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coupling between elements while the Mathcad application 
ignores these effects. 




FIGURE 7.3 Comparison of radiation patterns using ELNEC 
and Mathcad of a caged dipole in free space with L = A/2 
m, N = 8, and a = 0.125 m 
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B. CAGED DIPOLE ORIENTED VERTICALLY OVER EARTH 



The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented vertically 
over the earth. The height of the antenna above ground is a 
physical dimension obtained through photographs: 
h = antenna height above ground 
e r = relative dielectric constant of earth 
a = conductivity of earth 

av w = incoming wave electric field unit vector for 
vertical antenna 

0v p = antenna polarization direction for vertically 
orientated antenna 

As with the caged dipole in free space, the analysis of 
a caged dipole positioned vertically above the ground does not 
take into account the central conductor and support members 
since they produce negligible effects on the antenna's 
electrical properties. 

From image theory and pattern multiplication for arrays, 
the total far-field electric field intensity of a caged dipole 
oriented vertically above a flat earth is [Ref. 4: pp . 229- 

235] : 

£V c (0,<t>) =Ev(Q) AFv(6,<$>)(e jkhcosie) + r v e- jkhcos < e} ) ( V/m ) (7.28) 

In (7.28), Ev(0) is the far-field electric field intensity of 
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a finite length dipole in free space which is given by (7.5), 
and AFv (0,0) is the array factor for a circular array which is 
(7.1) . The total electric field intensity is valid only above 
the ground (0 < 0 < n/2) [Ref. 3: pp. 135-142]. The distance 
from the earth to the center of the antenna is (h) , and the 
vertical reflection coefficient (T v ) of the ground is [Ref. 4: 



pp. 231- 233]: 



rvo) 



e' r cos (0) - y/e ; r - sin 2 (0) 
t' z cos (6) + ^e^-sin 2 (0) 



( dimensionless ) (7.29) 



where 

e 7 = e - j — ( dimensionless ) (7.3 0) 

r z 2n fe 0 

In (7.29) and (7.30), (e' r ) is the relative complex effective 

dielectric constant of the ground, (e r ) is the relative 
dielectric constant of the ground, (a) is the conductivity of 
the ground, and e c is the permittivity of free space ( 1/(367 t) 
x 10" 9 ) . For a perfect ground plane, a -* °o and r v = 1 . 

The radiation intensity (U v (0,0)) of a caged dipole over 
earth is [Ref. 3: p. 28]: 

r 2 

U V (Q, <t>) = — — [|£v c (0, <J>) | 2 ] ( W/ solid angle) (7.31) 

The radiated power (Pv rad ) is present only above the earth. 
Thus, [Ref. 3 : p. 28] : 
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(7.32) 



Pv rad = f 2 *[>uv(d,<t>)sin(d)ddd<t> (W) 

Jo Jo 

Directivity (Dv 0 ) and effective isotropic radiated power 
(EIRPv) for a vertically oriented caged dipole over the earth 
are given by (7.12) and (7.13), respectively, with the results 
of (7.32) inserted. Uv max in (7.12) is simply the maximum 
value of (7.31) which, in Mathcad, is found by constructing 
(7.31) as a vector and then finding the maximum value of the 
vector . 

Input resistance (Rv in ) , voltage reflection coefficient 
(TV) , gain (Gh) , maximum effective aperture (Av em ) , and maximum 
effective height (hv em ) for a vertical caged dipole over the 
earth are determined using the equations for those parameters 
for a caged dipole in free space. 

The Polarization Loss Factor (PLFv) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (0v y ) is estimated. The antenna 
polarization vector (av a ) is then calculated and, estimating 
the unit polarization vector of the incoming wave (crv w ) , PLF 
is found from (7.21). 

Table 7.2 is a comparison of the gain of a vertical caged 
dipole over a flat earth obtained with both ELNEC and the 
Mathcad application. As can be seen, the gains are similar. 
As previously mentioned, the difference in the gains can be 
attributed to the fact that ELNEC takes into account mutual 
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coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 
and approximates ground effects. 



TABLE 7.2 Caged Dipole Antenna Gain Comparison 




Since the caged dipole behaves in a manner similar to 
that of a thin dipole except for bandwidth characteristics, 
the radiation pattern of a vertically oriented caged dipole 
over ground will resemble that of a thin vertical dipole over 
ground. Figure 7.4 is a plot of the radiation pattern of a 
caged dipole over earth. There is a very close resemblance to 
the radiation pattern of a thin vertical dipole over ground 
[Ref. 4: p. 234]. Figure 7.5 is a comparison of radiation 
patterns computed by ELNEC and Mathcad software. As expected, 
the patterns are similar. 
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FIGURE 7.4 E-plane radiation pattern of a caged dipole 
with L = A/2 m, N = 8, a = 0.125 m, h = 0.25 m, a = 0.01 

S/m, and e r = 15 




— Mothcod Pattern, ELNEC Pattern 



FIGURE 7.5 E-plane radiation pattern comparison between 
ELNEC and Mathcad of a caged dipole with L = A/2 m, N = 
8, a = 0.125 m, h = 0.25 m, a = 0.01 S/m, and e r = 15 
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C. CAGED DIPOLE ORIENTED HORIZONTALLY OVER EARTH 

The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented horizontally 
over the earth: 

h = antenna height above ground 

e r = relative dielectric constant of earth 

a = conductivity of earth 

ah w = incoming wave electric field unit vector for 
vertical antenna 

0h p = antenna polarization direction for vertically 
orientated antenna 

As with the evaluations of a caged dipole both in free 
space and oriented vertically over the earth, the analysis of 
a caged dipole positioned horizontally above the ground does 
not take into account the central conductor and support 
members since they produce negligible effects on the antenna's 
electrical properties. 

Through the use of image theory, the far-field electric 
field intensity of a caged dipole oriented horizontally above 
a flat earth can be obtained. The 0 component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 
space is [Ref. 3: pp. 143-144]: 
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( V/m) 



(7.33) 



£^(0,40 =jr i In& 






2 nr 



it 



cos|-^sin (0) sin(4>)j -cos|-^ 
\/l - sin 2 (0) sin 2 (4>) 



The 0 component of the total electric field over ground 
is [Ref. 4: pp. 229-235]: 

£0h t (0,<t>) =£0h(0,<t>) AFh(B,$){e jkhcosm - r v e^ khcos(6) ) {V/m) (7.34) 



The array factor (AFh(0,0)) in (7.34) for a horizontal array 
is determined by [Ref. 3: pp. 274 and 776]: 

AFh(B,<\>) = I n e jkasinesin * sin * n ( dimensionless ) (7.35) 

71 * 1 



In (7.33) through (7.35), the electric field (E0h(0,0)), the 
total electric field (E0h t (0,0)), and array factor (AFh(0,0)) 
are valid only above the ground (0 < 0 < n/2) [Ref. 3: pp. 
135-142]. In (7.34), (h) is the distance from the earth to 
the center of the antenna, and the vertical reflection 
coefficient (F v ) of the ground is given by (7.29). 

The $ component of the electric field of a horizontal 
dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 
135-142] : 



F4>h(0,<j)) =JTl 0 ^ e 






2nz ff 



cos|-^sin (0) sirx4>)|- cos|— 
v/l — sin 2 (0) sin 2 (({>) 



(V/m) 



(7.36) 



From image theory, the $ component of the total electric field 
over earth is [Ref. 4: pp. 229-235]: 
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E<\>h c (6,4)) =Ebh(B,4>) AFh(Q, 4>)(e^ cos(0 > 



+ r h e' jkhcosl6) ) ( v/m) (7.37) 



In (7.37), (AFh(0,0)) is given by (7.35), and (T h ) is the 

horizontal reflection coefficient of the ground and is [Ref. 
4: pp. 231-232]: 



cos(0) - y/e^- sin 2 ( 0 ) 
cos(0) + yjt r r - sin 2 ( 0 ) 



( dimensionless ) (7.38) 



When a -* «, r h = -1. 

The radiation intensity (U h (0,</>)) of a caged dipole over 
a earth must take into account both the 0 and $ components of 
the electric field. Therefore, the radiation intensity is 
[Ref. 3: p. 28]: 

2 

Uh(6,<p) = — — [|£0h t (0, 4>) \ 2 + \E$h c (0, (J>) | 2 ] ( w/ solid angle) (7.39) 

2 , n 0 

The radiated power (Ph rad ) is present only above the earth. 
Thus , 

p hrad = f f ' Uh ( 0, 4> ) sin ( 0 ) dQ c?4> (W) (7.4 0) 

Jo J o 



The remaining parameters (i.e., Dh 0 , EIRPh, Rh in , etc.) 
for a horizontal caged dipole over the earth are calculated 
using the equations for these parameters listed previously. 

Table 7.3 is a comparison of the gain of a horizontal 
caged dipole over a flat earth obtained with ELNEC and the 
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Mathcad application. As can be seen, the gains are similar. 
As previously mentioned, the difference in gains can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 
and approximates ground effects. 



TABLE 7.3 Caged Dipole Antenna Gain Comparison 



ANTENNA 

DIMENSIONS 


ELNEC 
Gain (G) 


Mathcad 
Gain (G) 


L - A/2 m 






a — 0.025 m 




| 


N = 8 
h = 1.0m 
a = 0.01 S/m 
e r = 15 


4.99 dB 

L 


4.73 dB 



As previously mentioned with regard to the vertically 
oriented caged dipole over the earth, the caged dipole behaves 
in a manner similar to that of a thin cylindrical dipole. 
Therefore, the radiation patterns of a horizontally oriented 
caged dipole are similar to those of horizontal thin 
cylindrical dipole which in turn resembles that of a thin 
dipole. Figures 7.6 and 7.7 illustrate the H-plane and E- 
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plane patterns, respectively, of a caged dipole over the earth 
[Ref. 4: pp. 231-237], The radiation patterns of the 
horizontally orientated caged dipole over a flat earth are 
very similar to those of a horizontal thin dipole over earth. 
Figures 7.8 and 7.9 compare the radiation patterns between 
those computed with ELNEC and those computed by Mathcad 
software. As expected, the patterns are similar. 




FIGURE 7.6 H-plane radiation pattern of a caged dipole 
with L = X/2 m, N = 8, a = 0.125 m, h = 1.0 m, a = 0.01 
S/m, and e r = 15 
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FIGURE 7.7 E-plane radiation pattern of a caged dipole 
with L = X/2 m, N = 8, a = 0.125 m, h = 1.0 m, a = 0.01 
S/m, and e r = 15 




FIGURE 7.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a caged dipole with L = X/2 m, N = 
8, a = 0.125 m, h = 1.0 m, a = 0.01 S/m, and e r = 15 
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FIGURE 7.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a caged dipole with L = X/2 m, N = 
8, a = 0.125 m, h = 1.0 m, a = 0.01 S/m, and e r = 15 



84 



VIII. PARABOLIC REFLECTORS 



Many communications and radar operations require very 
large values of antenna gain (G) that are difficult to achieve 
with a single device. One of the most popular ways to 
construct a very high gain antenna is to reflect the 
electromagnetic energy of a relatively small antenna off a 
metallic parabolic dish in the far-field of the feed device 
[Ref. 13: p. 12-2]. If the source antenna is located at the 
focus of the parabola and if the radius (a) of the parabolic 
dish periphery is large with respect to wavelength (A) , then 
the reflector may be approximated as a uniform phase aperture 
antenna. Gain in excess of 30 dB is common for this type of 
structure [Ref. 4: pp. 423-424]. The geometry of a parabolic 
reflector is illustrated in Figure 8.1, where the parameters 
include the focal distance (f), the distance (p) from the 
origin to any point (P) on the reflector surface, and the 
distance (r') from the origin to the projection of point (P) 
onto the z = 0 plane. The primed angles in Figure 8.1 
correspond to the parameters associated with the feed horn. 
The Mathcad application user should note that the focal 
distance should exceed the minimum distance to far-field for 
the feed antenna. 
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The Mathcad parabolic reflector applications are unique 
in that the software is not written as a stand alone package. 
Rather, the reflector software is run as an addendum to the 
programs that already exist for various antennas [Ref. 2: pp. 
4-15, 67-114]. Four types of feeds are analyzed: helical, 
planar equiangular spiral, pyramidal horn, and conical horn. 
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The following additional inputs are required to analyze 
parabolic reflectors: 

foe = reflector focal length 

a ref = radius of the mouth of the reflector 

E 0 = electric field scale factor 

e b = blockage efficiency 

e sp = spar efficiency 

r ff2 = secondary far field observation distance 
x,y,z = polarization loss factor coordinates 
o w = incoming wave unit polarization vector 
tl, t2 = radiated power increments for (p and 6 
il = secondary field increments 
N , M = summation increments 

The first two inputs are physical dimensions obtained 
primarily through photographs, inputs with the * are 
parameters that are either known or estimated, and the last 
inputs are only used to affect the processing time of the 
application . 

In order to calculate the far-field patterns of a 
parabolic reflector, the electric field at the aperture of the 
device must be known. Since the electric field patterns of 
interest are in the far-field, it is assumed that the circular 
aperture created by the reflector is located at the origin and 
is parallel to the z = 0 plane. The Cartesian components of 
the aperture electric field are expressed as [Ref. 14: pp. 
121-123; Ref. 15: pp. 29-31]: 
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(V/Ttl) 



( 8 . 1 ) 



E ax (r 7 ,0 7 ) 



[- E ef (0 7 (r 7 ) ,0 7 )cos(0 7 ) 

+ E^ f (0 7 (r 7 ) ,0 7 ) sin(0 7 ) ]e 



jk (p(r' ) -2foc ) 



E ay (r 7 ,0 7 ) = [ E ef (0 7 (r 7 ) ,0 7 ) sin(0 7 ) 



+ E^ £ (0 7 ( r 7 ) , 0 7 ) cos (0 7 ) ] e j k(p(r,) ' 2foc) (V/m) 



( 8 . 2 ) 



In (8.1) and (8.2), (E ax (r',0')) and (E ay ( r' , 0' ) ) are the 

aperture electric fields and k is the wavenumber (2ir/X). 

The relationships between (0' ) , (p) , and (r 7 ) needed to 

evaluate (8.1) and (8.2) are expressed as [Ref. 4: p. 426]: 



0 7 = 2tan 1 



2 foe 



{radians) 



(8.3) 



o(r') = 4 ( foe) 2 + (.r 7 )^ 
4 foe 



(w) 



(8.4) 



The electric field from the entire parabolic assembly is 
commonly referred to as the secondary pattern. Referring to 
Figure 8.1, the electric field in the -z half-plane is 
considered the secondary pattern. The secondary field of a 
parabolic reflector is calculated using the electric field 
integral equation solution and computing the component vector 
potentials (P x (0,0)), (P y (0,0)) as follows [Ref. 4: p. 426]: 



w -rr 



E (r / ,<{> / ) [e-?** r/sin ( 0 )(cos ($/) cos (<t>) 



sin (<j>') sin (<J>) )j 



*r'dr , d<b' 



( V/m) 



(8.5) 
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( 8 . 6 ) 



2?r na rp f 

p 5( 0 .« ) =J o | o E„(r',0')[e 




* r / dr ' d0 f 



(V/m) 



Double integrals in Mathcad take significantly longer to 
converge to an answer than do single integrals [Ref. 15: p. 
200]. To reduce the processing time of the applications, 
(8.5) and (8.6) were each separated into two single integrals. 
Therefore, for the applications, (8.5) and (8.6) are each 
computed as: 



A (0 , <J) , 4>0 =f a '* f E a (r 4)0 e J>cr ' sin,e) t cos( ' t,,) c os ' , t >) * sinful sin (4»] r / dr i (8.8) 



Any contribution from feed system back lobes is ignored, 
and the component vector potentials calculated in (8.7) and 
(8.8) are applied to the subsequent formulas to compute the 
secondary field [Ref. 4: p. 383]: 




( V/m) 



(8.7) 



where 




[P x (0, 4)) cos (4>) + P y (0, 4>) sin (<t>) ] (V/m) (8.9) 




[-P x (0,<|)) sin(<J>) + P y (0, <{>) cos (<t>) ]cos (0) (8.10) 
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The radiation intensity (U(0,0)), radiated power (P rad ) , 
and directivity (D 0 ) of the secondary field pattern are 
expressed using the intrinsic impedance of free space ( r? 0 = 
1207T) as follows [Ref. 3: pp. 27-30]: 

t/(0,4>) = [[£■„ (6, <t>) | 2 + 1-Si (0, 4>) | 2 1 (W/ solid angle) (8.11) 

2ri 0 1 

P rad = f o 2n f o n/Z U(6,(t>)sin(d)ddd < p (W) (8.12) 

D _ 4 71 (dimensionless) (8.13) 

° P 

rad 

In (8.13), (U raax ) is the maximum value of radiation intensity, 
obtained by creating (U max ) as a matrix and finding the maximum 
value of the matrix. 

The computer processing time required to evaluate (8.12) 
is prohibitively long. Thus, the parabolic reflector Mathcad 
applications use the following numerical approximation to 
evaluate radiated power [Ref. 3: pp. 37-42]: 




In (8.14), (N and M) are user selected integers. This 
approximation for P rad reduces the processing time of the 
applications substantially. As N and M increase, the results 
obtained with (8.14) converge to that obtained with (8.12). 
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When N = M = 3, the processing time for only P rad of the 
helical feed reflector option is approximately four hours on 
a 33 MHz, 386 machine, and when N and M is increased to 5, the 
time to calculate P rad for the same application is ten hours. 
For the applications, we chose N = M = 10 to evaluate P rad with 
reasonable accuracy. With this choice, it takes approximately 
one day for the application to run, but this is much less than 
the thirty-six hours it takes if (8.12) is used. 

Parabolic reflectors are members of the aperture antenna 
family. The gain of an aperture antenna is the product of its 
aperture efficiency (e ap ) and directivity. Therefore, the gain 
of a parabolic reflector is written as [Ref. 6: p. 575]: 

G = t ap D 0 ( dimensionless ) (8.15) 

Aperture efficiency is the product of several separate 
terms, including spillover efficiency (e s ), taper efficiency 
(e t ), phase (or random surface error) efficiency (e p ), 
polarization efficiency (e x ) , blockage efficiency (e b ) , spar 
efficiency (e sp ) , and ohmic efficiency (? 0 hmic) • 

Spillover efficiency (e s ) is the fraction of power 
radiated by the feed system that is intercepted and reflected 
by the parabolic dish. Spillover efficiency is the most 
important efficiency term for any feed system and is defined 
as [Ref. 6: pp. 583]: 
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( dimensionless ) (8.16) 



f 2 " f° 0 U f (Q') sin (00 dO'dc}/ 
_ J o Jo 

rfuAO') sin (00 dQ' d$' 

Jo Jo 



In (8.16), (U f (0')) is the radiation intensity function of the 
feed antenna and (0 O ) is the value of (8.3) when (r' ) equals 
the radius of the mouth of the reflector (a ref ) . 

Taper efficiency (e t ) accounts for the lack of amplitude 
uniformity of the feed pattern on the surface of the 
reflector. Taper efficiency, which is a dimensionless 
quantity, is expressed as [Ref. 3: pp. 626-627]: 



e 



= 2cot 2 





L 


r 2n 

Jo 


r° 

Jo \ 


4nU f (d') | 

L I - ^ n 


2 

(IH*' 


U Cli 1 

Prad 1 


2 J 


/, 


2 * [8° — sin (00 dQ'dfy' 

3 Jo P iad 



(8.17) 



Phase efficiency (e p ) is a function of the power loss 
that occurs if the field at the mouth of the dish is not in 
phase at every point in the aperture. Phase efficiency is 
approximately given by [Ref. 4: pp. 434-435]: 

-47i6-10“ 5 a.. / r 

I — , ,. 7 , (8.18) 

e p = e A (dimensionless) 

The power that is lost due to cross-polarized electric 
fields in the aperture of the reflector determines 
polarization efficiency (e x ). Polarization efficiency is 
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difficult to precisely calculate, but can be reasonably 
approximated as 0.98 [Ref, 4: p. 435]. 

Blockage (e b ) and spar (e sp ) efficiencies are determined 
by the presence of feed structures, support struts, and 
possible signal processing eguipment in front of the 
reflector's aperture. Precise computation of blockage and 
spar efficiencies requires an extensive method of moments 
simulation and is not attempted by the Mathcad parabolic 
reflector applications. However, if the radius (a f ) of the 
blocking structure at the focal point and the number of spars 
(N) that are (1/2) thick is known, blockage efficiency and 
spar efficiency can be estimated using Tables 10.1 and 10.2, 
respectively [Ref. 4: p.436]. 



TABLE 8.1 Blockage Efficiency 





0.05 


0 . 10 


0.20 




0.99 


0.96 


0.84 



TABLE 8.2 Spar Efficiency 





a re f = 5/1 


a ref = 50/1 


a re£ = 100/1 


N = 3 


0.95 


0.99 


1.00 


N = 4 


0.94 


0.99 


1 . 00 
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Ohmic efficiency (e ohmic ) for a reflector is usually very 
high but difficult to accurately calculate. The Mathcad 
parabolic reflector applications assume ohmic efficiency is 
0.98. 

With all reflector efficiency terms known or estimated, 
the aperture efficiency (e ap ) , effective aperture (A e ) , and 
isotropic radiated power (EIRP) are given by [Ref. 3: pp. 623- 
630] : 

e ap = e s e t e p e * e * c sp c ohmic ( dimensionless ) (8.19) 

CapM^W) 2 (/n2) (8.20) 

EIRP = P zad D 0 ( W ) (8.21) 



The magnitude of the feed current (I Q ) for the source 
antenna cannot be estimated based on measured geometry. 
However, if the feed current is normalized to 1 amp, radiation 
resistance (R r ) for the entire reflector structure is: 



R 



r 



2 £Vad 

Uol 2 



(ft) 



( 8 . 22 ) 



In (8.22), (P rad ) is the radiated power from the aperture of 
the reflector, not of the feed antenna. 

The Mathcad parabolic reflector applications use the 
radiation resistance computed for the parabolic reflector 
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antenna system to determine effective height (h em ) by [Ref. 6: 
p. 42 ] : 




(8.23) 



The bandwidth of a parabolic reflector is assumed to be 
the same as its feed antenna. 

The polarization loss factor (PLF) of the parabolic 
reflector for a user selected point in the secondary far-field 
and a specific wave unit polarization vector (a w ) is found 
using the electric field components of (8.9) and (8.10) as 
follows [Ref. 3: 51-53]: 



and the Cartesian components of the electric field are [Ref. 
3: p. 776]: 



PLF = |o w .-o*| 2 (dimensionless) 



(8.24) 



where the unit polarization vector of the antenna is: 



a a (x,y, z) = 



^_x_ + ^ y + ^_z 



( dimensionless ) (8.25) 



v/j E(x,y, z) I 



E x = E e cos (0) cos (4>) - E 4> sin(<})) (v/m) 



(8.26) 



E y = E e cos (0) sin (4>) + E^c os(4>) (V/m) 



(8.27) 



E z = - E e sin (0) ( v/m) 



(8.28) 
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All of the secondary field parameters calculated by the 
Mathcad parabolic reflector applications assume that the 
observation point is in the far-field. Thus, all of the 
subsequent formulas must hold for the application calculations 
to be valid [Ref. 4: pp. 24-25]: 

r ^ 1.61 (m) (8.29) 



r > 10 a ref (m) 



(8.30) 



r 2 i 



8( a ref) 2 

X 



(dimensionless) 



(8.31) 



Table 8.3 is a comparison of reflector directivity 
computed by the Mathcad applications to predictions made for 
an identical antenna with unity aperture efficiency [Ref. 3: 
p. 634]. As can be seen, the comparisons vary by 
approximately 8 dB. This can be attributed to the fact that 
P rad was approximated using (8.14) instead of (8.12), and also 
to the fact that the numerical tolerance variable (TOL in 
Mathcad's math menu) for the application was increased to 0.5 
verses 0.001 in order to reduce processing time. This 
increase in TOL reduces both the accuracy of (8.14) and the 
processing time of the application. 
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TABLE 8.3 Parabolic Reflector Data Comparison 



ANTENNA 


PREDICTED 


CALCULATED 


PARAMETER 


DATA 


DATA 

i 


DIRECTIVITY 


42 dB 


50.2 dB 



Figure 8.2 is a comparison of measured data fed by an 
optimum horn with that calculated by the Mathcad application 
fed by an equiangular spiral for a parabolic reflector (a ref = 
0.320 m, f = 0.147 m) operating at 35 GHz [Ref. 14: p. 126]. 
As can be seen, the comparison is good. The Mathcad pattern 
of Figure 8.2 is not symmetric which can be attributed to the 
number of increments used to increase the speed of processing. 

The reflectors fed by pyramidal and conical horn antennas 
are not included in this report. The applications have been 
developed and are included on the diskette containing the 
software. The complexity of the software results in 
prohibitively long processing time for the two applications 
(i.e., the pyramidal horn fed reflector application was still 
calculating P rad2 after seven days when run on a 50 MHz, 486 
personal computer using Mathcad 4.0). 
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FIGURE 8.2 Comparison of Parabolic Reflector E-Plane 
Electric Field Patterns 
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IX. REMARKS AND CONCLUSIONS 



This thesis and the accompanying software applications 
are intended to provide the Naval Maritime Intelligence Center 
(NAVMARINTCEN) with software for the analysis of various 
antenna types where information regarding the antenna is 
restricted primarily to the antenna type and the physical 
dimensions of the antenna. The goal of this research was to 
provide user friendly software that an engineer familiar with 
basic antenna types can easily use and interpret. Although 
the software was developed for that purpose, the complexity of 
formulas that apply to many of the antenna types of necessity 
reduced the simplicity of the applications. 

The software applications developed are written in 
Mathcad 3.1 and are compatible with any IBM personal computer 
that supports Mathcad 3 . 1 or Mathcad 4.0 for Windows. Many of 
the applications require extensive processing time. In order 
to reduce processing time, which results from the complexity 
of various integrals, approximate numerical techniques were 
implemented to evaluate the integrals and the numerical 
tolerance of the more complicated applications was reduced. 
Even with the processing time reduced, some programs require 
up to five days to complete when run with Mathcad 4.0 on a 50 
MHz, 486 personal computer. Applications that did not exceed 
a day to complete were run on a 33 MHz, 386 personal computer 
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with Mathcad 4.0. All applications were run using Mathcad 4.0 
since Mathcad 4.0 is approximately twice as fast as Mathcad 
3.1. As previously mentioned, the applications were written 
in Mathcad 3.1 due to sponsor requirements. Mathcad 3.1 is 
compatible with Mathcad 4.0, but Mathcad 4.0 is not compatible 
with Mathcad 3.1. 

For this thesis, a number of the antennas could not be 
adequately researched through the use of texts. In some 
cases, such as the caged dipole, very little information on 
the antenna could be found in either texts or in professional 
journals. For the caged dipole, the antenna is modeled as a 
circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles to take into account the number of 
conductors that are present in the caged dipole. This model 
provides good results when compared to predicted (the caged 
dipole's performance is similar to that of a thick cylindrical 
dipole) and ELNEC results. 

It was found that Mathcad is not the ideal tool to 
analyze antennas due to its long processing time. An 
additional drawback of Mathcad is that when a function is in 
the beginning of an application and the same function is 
placed later in the application, Mathcad does not store the 
values of the function, but re-evaluates the function. This 
re-evaluation increases processing time significantly, 
especially when the function itself takes a very long time 
(hours vice minutes) to reach a numerical solution. 
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Nevertheless, the Mathcad applications developed in this 
research allow the user to analyze several antenna types 
(i.e., arrays with isotropic point sources, parabolic 
reflectors with the various feed antennas) that could not 
easily be evaluated using some antenna numerical analysis 
programs such as NEC, ELNEC, or MININEC. For the other 
antenna types researched, the antenna numerical analysis 
programs mentioned would be easier to setup, run, and provide 
parameters closer to those of a real antenna. The 
disadvantage of this approach is that it requires an ability 
to evaluate antennas using a program such as ELNEC, while the 
Mathcad applications do not require any advanced programming 
knowledge on the part of the user. 
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LINEAR ARRAY 

MATHCAD SOFTWARE ARRAY_LN . MCD 

This linear array application is based on N isotropic point sources 
uniformly spaced and aligned on the z-axis. Mutual interference between 
adjacent sources are ignored when calculating antenna parameters. Since 
isotropic point sources are used as the radiating elements, polarization 
of the antenna can not be determined. The gain can only be idealized 
when it is assumed that the antenna input resistance is equivalent to 
the radiation resistance. Bandwith calculations are not carried out as 
a result of the narrowband characteristics of arrays, and in large part 
to the dependance inter-element spacing (d) has with the wavelength (1) . 
The following limits should be adhered to when analyzing Broadside, 
End-Fire, and Phased (scanning) antennas: 

0 < 0 o < 7 i (Scanning Antennas) 

0 O = 7i/2 (Broadside Antennas) 

G 0 = 0 or n (Ordinary End-Fire Antennas) 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn ) . ) 

The linear array antenna Mathcad application will compute the following 



parameters (Items with * indicate parameters 


that can be computed but 


their 


results are idealized as the result of 


using isotropic point 


sources as the radiating elements) : 




?. = 


Wavelength 




r min 


= Minimum Distance to the Far-Field 




P = 


Phase Difference Between Adjacent Elements 


V = 


Array Factor Phase Shift 




AF n 


= Normalized Array Factor 




U(0) 


= Radiation Intensity 




p rad 


= Radiated Power 




Do = 


Di rectivi ty 




EIRP 


= Effective Isotropic Radiated Power 




R r = 


Radiation Resistance 




R in 


= Antenna Input Resistance * 




r = 


Voltage Reflection Coefficient * 




£ t = 


Antenna Efficiency * 




G = 


Gain * 






= Maximum Effective Aperture * 




^em 


= Maximum Effective Height * 
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The following known or estimated data must be entered: 



N = Number of Isotropic Radiating Elements 

f = Frequency of Interest 

d = Distance Between Adjacent Elements 

0 O = Direction of Main Lobe 

I Q = Antenna Feed Current 

Z Q = Characteristic Feed Impedence 



N := 10 



Enter input data here 



( elements ) 



d :=0 25 



f : = 3 1 0 8 (Hz) 



0o := 0 



Io := 10 



(A) 



Zo-15.0 



Define constants and calculate wavelength: 



c- 2.9979 10 8 (meters/sec) 



i]o •= i 20 7i 



c 

f 



^ 71 

(meters/cycle ) k := — 

). 



1=0.999 



(meters/cycle ) 



ccd •= 1 



(meters ) 



( radians ) 



(Q) 



(Q) 



(m 1 ) 



(dimensionless ) 
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Calculate linear array parameters: 
Define angular offset 0: 



0 := (o- ] 0 ‘ 4 ) , (— ) . 

V 1 80/ 


(2 71- 10‘ 4 ) 


( radians ) 


Minimum Distance 


to Far-Field r mln : 




V s1 - 6 '* 


(meters ) 




r, :=5 N d 


(meters ) 




2 (N d) 2 

V = . 


(meters ) 




rmin := max(r) 


(meters ) 




rmin = 12.509 


(meters ) 




Array Factor AF n 







P - k d cos(0o) (radians) 



y(0) := k d cos(0) -t p 



( radians ) 



AFn(O) 



N 



sin| i|/(0) 
sin|i- iji(0) 



(dimensionless ) 
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Radiation Intensity U: 



U(O): = (AFn(0)) 2 



(W/solid angle) 



Radiated Power P r ^ H : 



Prad '= 2-7i- 



U(0) sin(G) d() 



(W) 



Prad = 1.256 



(W) 



Directivity D n : 



Umax 1 



(W /solid angle) 



Do := 



4 7i Umax 



Prad 



(dimensionless ) 



Do = 10.007 



( dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP := Prad Do 



(W) 



EIRP = 12.566 



(W) 
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Radiation Resistance R r : 



2 Prad 

C|lo| ) 2 



(Q) 



Rr = 2.5 1 2 



(Q) 



Input Resistance R in : 



Rin := Rr 



(Q) 



Rin = 2.5 1 2 



(Q) 



Voltage Rejection Coeffecient F: 



Rin - Zo 
Rin + Zo 



(dimensionless ) 



f =“0.935 



(dimensionless ) 



Reflection Efficiency : 

et = ecd [ 1 - ( | r| ) 2 ] ( dimensionless ) 



st =0.125 



(dimensionless ) 
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Gain G: 



G- = et Do (dimensionless) GdB •= 10-log(G) 



G = 1.255 (dimensionless) GdB =0.986 



Maximum Effective Aperture 



Aem := 



Gr 

4*71 



Aem = 0. 1 



Maximum Effective Height h Arn : 



hem := 2 



r RrAern 



T]0 



hem = 0.052 



(dB) 



(dB) 



(m 2 ) 



(m 2 ) 



(m) 



(m) 
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NORMALIZED ARRAY FACTOR RECTANGULAR PLOT 



Number of radiating elements: 



N = 10 



Direction of Main Lobe (radians): 

Go = 0 




0 0.628 1.257 1.885 2.51? 3.142 3.77 4 398 5.027 5.655 6.28? 

0 
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NORMALIZED ARRAY FACTOR POLAR PLOT 



Number of Elements: 
Direction of Main Beam 



( radians ) : 



N = 10 
Oo = 0 



X(0) := | AFn(0)| 



COS ( 0+ ^J 



Y(0) 




-1 



0 = 0 ° 




0 = 180 ° 
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PLANAR ARRAY 

MATHCAD SOFTWARE ARRAY_PL . MOD 

This planar array application is based on M and N isotropic point 
sources uniformly spaced and aligned on the x-axis and y-axis 
respectively. Mutual interference between adjacent elements are 
neglected when calculating antenna parameters. Since isotropic point 
sources are used as the radiating elements, polarization of the antenna 
can not be determined. The gain can only be idealized when it is 
assumed that the antenna input resistance is equivalent to the radiation 
resistance and the total antenna efficiency is calculated for a ideal 
lossless antenna. 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn) . ) 

The planar array antenna Mathcad application will compute the following 



parameters (Items with * indicate parameters 


that can be 


computed but 


their 


results are idealized as the result of 


using isotr 


opi 


c point 


sources as the radiating elements) : 










y — 


Wavelength 










L x = 


Antenna Dimension in x-direction 










L y = 


Antenna Dimension in y-direction 










r min 


= Minimum Distance to the Far-Field 










Px = 


Progressive Phase Shift Between Adjacent 


Elements 


in 


x-direction 


"CD 

II 


Progressive Phase Shift Between Adjacent 


Elements 


in 


y-di rection 


Yx = 


Array Factor Phase Shift in x-direction 








Yy = 


Array Factor Phase Shift in x-direction 








AF n (6, 


<j)) = Normalized Array Factor 










U 


= Radiation Intensity 










p rad 


= Radiated Power 










Do = 


Directivity 










EIRP 


= Effective Isotropic Radiated Power 










R r = 


Radiation Resistance 










R m : 


= Antenna Input Resistance * 










r = 


Voltage Reflection Coefficient * 










£ t = 


Antenna Efficiency * 










G = 


Gain * 












= Maximum Effective Aperture * 










^em 


= Maximum Effective Height * 
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The following known or estimated data must be entered: 



M = Number of Isotropic Radiating Elements in x-direction 
N = Number of Isotropic Radiating Elements in y-direction 
f = Frequency of Interest 

d x = Distance Between Adjacent Elements in x-direction 

dy = Distance Between Adjacent Elements in y-direction 

0 O = Direction of Main Beam at 0 = Oo 

(j) 0 = Direction of Main Beam at = 4>o 

I Q = Antenna Feed Current 

Z Q = Characteristic Feed Impedence 

j = Number of Increments in Degrees for Far-Field Radiation Pattern 



Enter input data here 



M := 5 



(elements ) 



d\ := 0.5 



(meters ) 



N - 5 



(elements ) 



dv =0.5 



(meters ) 




(Hz) 



0o := 



n 



{ radians ) 



6 



Io =1.0 



(A) 



( radians ) 



4 



Zo i = 75.0 



(Q> 



180 



( degrees ) 
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Define constants and calculate wavelength: 



c := 2.9979 10* (meters/sec) 



t }o := 1 20 7i 



X:=‘ 


(meters/cycle ) 


k := 


f 






X =0.999 


(meters/cycle ) 


CCd 


Calculate 


linear array parameters: 




Define angular offset Gandjk 




0 : (0- 10‘ 4 ) 


,(- lff 4 +— 1 ( 271 - Iff 4 ) 






l j I 





2_n 

X 



4> : = (o — iff 4 ) io‘ J + — j (21- icr 4 ) 



Minimum Distance to Far-Field r 



mm_L 



i’ 0 : = 1.6 /. (meters) 



5 ^(M dx) 2 1 (Ndv) 2 



(meters ) 



r 3 : = 



_ 2 [^/(M dx) 2 + 



(N-dy)" 



( me t e r s ) 



nmin := max(r) (meters) 



nmn = 25.01 8 



(Q) 



(m 1 ) 



( dimensionless ) 



( radians ) 



( radians ) 



(meters) 
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Array Factor AF n : 

Px := - k- dx- sin( 0o) * cos( 4>o) (radians) 

Py *=-k-dy sin(0o) sin(<K>) (radians) 



iyx(0,<t>) *=k dx sm(0) cos(<t») -h Px 



( radians ) 



i|/y(0,<f>) ‘ = k dy sin(0) sin(<j») + p\ 



AFn(0,<(>) 



sin | 


^ ^rx(e, <t> )j 


r 


f / N 

sin — 
\2 




1 

1 


M sinH 0 '**] 


1 


N sin| 


l 2 ) 



Radiation Intensity U: 



U(0,<)>) .= AFn(0,<j>) AFn(0,(j)) 



( radians ) 



(dimensionless ) 



(W/solid angle) 



Radiated Power P 



rad_l 



Prad := 



U(0.<j>) sin(0) dO d4> 



(W) 



Prad =0.206 



(W) 
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Directivity D n : 



Umax := 1 



(W /solid angle) 



Do := 



4 7i- Umax 



Prad 



(dimensionless ) 



Do =61.09 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP = Prad Do 



(W) 



EIRP = 12.566 



<W) 



Radiation Resistance R r : 



Rr : = 



2 Prad 

(llol ) 2 



(fl) 



Rr = 0411 



(Q) 



Input Resistance R n n : 

Rin := Rr (Q) 



Rin =0.41 1 



(Q) 
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Voltage Relection Coeffecient Y: 



_ Rin - Zo 
Rin + Zo 


(dimensionless ) 


r =-0.989 


(dimensionless ) 



Reflection Efficiency z t : 



ct := ecd [ 1 


- (|r|) 2 ] 




(dimensionless) 


el =0.022 






(dimensionless ) 


Ga in G i 








G := ct Do 


( dimensionless ) 


GdB := K) lop(G) 


(dB) 


G = 1.326 


(dimensionless ) 


GdB = 1 225 


(dB) 


Maximum 


Effective Aperture : 






n , 2 

A G /. , O 

Aem ■= (m z ) 

4-ji 


Aem =(J 105 


(m^ ) 


Maximum 


Effective Height h^ m : 








— 






hem := 2 J 


*Rr Aem 

(in) 

TP 


hem = 0 02 1 


(m) 
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TWO-DIMENSIONAL ANTENNA PATTERN: (x-z) plane 



Number of 


Elements 


along x-axis: 


M =5 




Number of 


Elements 


along y-axis: 


N =5 




Direction 


of Main 


Beam: 


0o = O.524 


( radians ) 








o 

II 

o 

-~J 

oc 


( radians ) 



X1(0) '= |AFn(0,O)| cosfo-— j Y1(0) ’= | AFn(0,O)| sin^0 — 

0 = 0 ° 
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TWO-DIMENSIONAL ANTENNA PATTERN: (y-z) plane 



Number of Elements along 
Number of Elements along 
Direction of Main Beam: 



x-axis : 
y-axis : 



X2(0) 



AFn 




cos 0- 



M =5 




N =5 




0o=O.524 


( radians ) 


<j)o = 0.785 


( radians ) 



Y2(0) 



AFn 





0 = 0 ° 




xrel . \2( 0) 
m ' 
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TWO-DIMENSIONAL ANTENNA PATTERN: $ = n/ 4 



Number of Elements along 
Number of Elements along 
Direction of Main Beam: 



x-axis : 
y-axis : 



X3(0) 



AFn 




cos 




M =5 
N =5 



Oo =0.524 
4>o =0.785 



Y3(0) 




( radians ) 
( radians ) 




0 = 0 ° 
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TWO-DIMENSIONAL ANTENNA PATTERN: 



Number 


of 


elements 


along x-axis: 


2 

H 




Number 


of 


Elements 


along y-axis: 


N =5 




Direction 


of Main 


Beam: 


0o= 0.524 


{ radians ) 










O 

II 

GC 

Vi 


( radians ) 



X4(G) := | AFn(0, 4>o)| cosfo — • 



Y4(0) •= | AFn(0,itio)| sin 




0 = 0 ° 
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CIRCULAR ARRAY 

MATHCAD SOFTWARE ARRAY CI.MCD 



This circular array application is based on N isotropic point sources 
uniformly spaced on the x~y plane along a circular ring of constant 
radius. Mutual interference between adjacent elements are neglected 
when calculating antenna parameters. Since isotropic point sources are 
used as the radiating elements, polarization of the antenna can not be 
determined. The gain can only be idealized when it is assumed that the 
antenna input resistance is equivalent to the radiation resistance and 
the total antenna efficiency is calculated for a ideal lossless antenna. 



(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn ) .) 



The circular array antenna Mathcad application will compute thefollowing 
parameters (Items with * indicate parameters that can becomputed but 
their results are idealized as the result of usingisot ropic point 
sources as the radiating elements): 



/. = Wavelength 

r min = Minimum Distance to the Far-Field 
a n = Phase Excitation of the nth Element 
AF (0, <|>) = Normalized Array Factor 
U(0,<t>) = Radiation Intensity 

p rad = Radiated Power 
D Q = Directivity 

EIRP = Effective Isotropic Radiated Power 
R r = Radiation Resistance 

R in = Antenna Input Resistance * 

T = Voltage Reflection Coefficient * 
c t = Antenna Efficiency * 

G = Gain * 

A em = Maximum Effective Aperture * 
h em = Maximum Effective Height * 
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The following known or estimated data must be entered: 



N = Number of Isotropic Radiating Elements 
f = Frequency of Interest 
a = Radius of circle 

0 O = Direction of Main Beam at 0 = 0o 

<j> 0 = Direction of Main Beam at § - 

I Q = Antenna Feed Current 

Z Q = Cha racteris tic Feed Impedence 

j = Number of Increments in Degrees for Far-Field Radiation Pattern 



Enter input data here 



N := 10 



(elements ) 



a := ■ 



10 



(meters ) 



f:=3 10 8 



(Hz) 



0o := 0 



( radians ) 



Io := 10 



(A) 



4)0 := 0 



( radians ) 



Zo — 75.0 



(Q) 



j - 360 



( degrees ) 



Define constants and calculate wavelength: 
c := 2.9979 10 8 (meters/sec) r|o'=120;i 



L 



c 

f 



(meters/cycl e ) 




/. 



( 0 ) 



(m 1 ) 
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Define constants and calculate wavelength: 

). =0.999 (meters/cycle) red := I (dimensionless) 

Calculate linear array parameters: 

Define angular offset 6_ and _tj>: 



e:=( 


o- ior‘j) ,f- io- 4 + — 


(2 7T- Iff 4 ) 


( radians ) 




\ j 






♦ =( 


0- Iff 4 ) lff 4 + — 


I (2 a- Iff 4 ) 


( radians ) 




\ J 1 


r 





Minimum Distance to Far-Field r m in : 

i* 0 =1.6 ). (meters) r =5 (2 a) (meters) 



.. 2 (2 a) 2 

2 ' , 

). 



(meters ) 



rmin := max(r) (meters) 



nmn = 20.278 (meters) 



Array Factor AF: 

n = 1 .. N (elements) (]) (radians) 

n N 

a n := - k a sin(Oo) cos^o- (1>J| (radians) 
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Array Factor AF: 



AF(0,<j>) - = — 
N 



j a sin( 0) cos($- + ci^l 



( dimensionless ) 



Radiation Intensity U: 



U(0,*) := AF( 0 . 4>) AF( 0, <J>) 



(W/solid angle) 



Radiated Power P r ^d : 



2 x 



71 

i 



Prad := 



•0 



U(O,<t>) sin(0)d0 d<() 

• 0 



(W) 



Prad =0.536 



(W) 



Directivity D n ; 



Umax : = 1 



iW/solid angle) 



Do •= 



4 7i Umax 



Prad 



( dimensionless ) 



Do = 23.427 



{ dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP := Prad Do 



(W) 



EIRP = 12.566 



(W) 
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Radiation Resistance R r : 
2 Prad 



Rr : = - 



( |lo| ) 2 



(Q) 



Rr= 1.073 



(Q) 



Input Resistance R in : 

Rin := Rr (Q) 

Rin = 1.073 (Q) 



Voltage Relection Coeffecient £: 

^ Rin - Zo 
Rin + Zo 



( dimensionless ) 



r =-0.972 



(dimensionless ) 



Reflection Efficiency : 

et = ecd [ 1 - ( |r| ) 2 ] (dimensionless) 

et =0.056 (dimensionless) 
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Ga in G i 



G- = ctDo (dimensionless) GdB •= 101og(G) (dB) 



G = 1.303 (dimensionless) GdB = 1.149 (dB) 



Maximum Effective Aperture 

ri 2 

Aem := (m 2 ) Aem =0.104 (m 2 ) 

4-ti 



Maximum Effective Height h ^r 



hem- 2 | ^ cm ( m ) hem =0.034 (m) 

J T)0 
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TWO-DIMENSIONAL ANTENNA PATTERN: (x-z) plane 



Number of 


Elements in x-y plane: 


N = 10 




Direction 


of Main Beam: 


Go = 0 


( radians ) 






<(>o = 0 


( radians ) 



X 1 ( 6) := |AF(0,O)| cos|g- Y 1(G) := | AF(0,O)| sm|o- ^ 

0 = 0 ° 




xref , X 1 ( 0 ) 
m 
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TWO-DIMENSIONAL ANTENNA PATTERN: (y-z) plane 



Number of 
Direction 



Elements in x 
of Main Beam: 



y plane: 



X2(0) 




N = 10 
0o = 0 
4>o = 0 

Y2(0) 



( radians ) 
( radians ) 




0 = 0 ° 
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TWO-DIMENSIONAL ANTENNA PATTERN; = n / 4 



Number of Elements in x-y plane: 
Direction of Main Beam: 



X3(0) 




N = 10 
Oo = 0 
4 >o = o 

Y3(0) 



( radians ) 
( radians ) 




0 = 0 ° 




xnsf.X3(0) 

m 
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TWO-DIMENSIONAL ANTENNA PATTERN: $ = 



Number of 
Direction 



Elements m x-y plane: 
of Main Beam: 



X4(0) := |AF(0,*o)| cos 0- ^ 



N = 10 




Go = 0 


( radians ) 


4>o = 0 


( radians ) 



V4( 0) := | AF(0.<t*o)| sin 0- ^ 



0 = 0 ° 




xref . \4{0) 
m 
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FOLDED DIPOLE 

MATHCAD SOFTWARE FOLDED . MOD 

To provide impedance matching to coaxial input cables and maintain good 
directional pattern characteristics, a folded dipole is used instead of 
a single dipole element. The geometry of a folded dipole is such that 
the spacing between dipole elements (d) is much less than the wavelength 
(X) and the length of an element (L) . 

This Mathcad application calculates the parameters of a folded dipole 
aligned vertically on the z-axis and a folded dipole positioned 
horizontally above the earth. For the calculations, the antenna is 
assumed to behave as a cylindrical dipole with equivalent radius (a e ) . 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., Rj_ n will be written Rin).) 

The folded dipole antenna Mathcad application will compute the following 
parameters ( * indicates antenna parameters that are calculated for a 
folded dipole in free space and positioned horizontally over the earth) : 



). = Wavelength 

a e = Equivalent Radius 

R r = Radiation Resistance 

R^ = Input Resistance of Linear Dipole 

X r = Reactance of Linear Dipole 

= Input Reactance of Linear Dipole 
= Input Impedance of Linear Dipole 



z o = 


Characteristic 


Impedance o 


z t = 


Characteristic 


Impedance T 


z in = 


Characteristic 


Impedance 


z o = 


Maximum Current 


of Folded : 


r min 


= Minimum Distance to the 



E (0) = Electric Field * 

E n (0) = Normalized Electric Field * 

U(0) = Radiation Intensity * 

Umax = Maximum Radiation Intensity * 

P rac j = Radiated Power * 

D 0 = Directivity * 

EIRP = Effective Isotropic Radiated Power 
Rin = Antenna Input Resistance * 

f = Voltage Reflection Coefficient * 



Transmission 
f Folded Dipol 
Dipole 



Line 

e 
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£ t = Antenna Efficiency * 

o a = Antenna Polarization Vector * 

PLF = Polarization Loss Factor * 

G = Gain * 

A em = Maximum Effective Aperture * 
hem = Maximum Effective Height * 

e r p = Relative Complex Effective Dielectric Constant of 
T v = Vertical Reflection Coefficient of Ground Plane 
f h = Horizontal Reflection Coefficient of Ground Plane 



Ground 



The following known or estimated data must be entered: 



f 

a 

b 

d 

L 

h 

£ cd 

c r 

0 

r ff 

1 

°wh 

^ph 



Frequency of Interest 

Radius of Conductor 

Radius of Conductor 

Distance Between Conductors 

Length of Antenna 

Height of Antenna Above Earth 

Conduction/Dielectric Efficiency 

Relative Dielectic Constant of Ground 

Conductivity of Ground 
: Far Field Distance 

Number of Increments m Degrees for Far Field Radiation Pattern 
Incoming Wave Electric Field Unit Vector For Free Space Antenna 

INcoming Wave Electric Field Unit Vector For Horizontal Antenna 
Antenna Polarization Direction For Horizontal Antenna 
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Enter input data here 



1 

aw := 0 
,0, 


( dimensionless ) 


'O' 

cwh •= 0 

111 


(dimensionless ) 


00 

o 

cb 

ii 


(Hz) 


Oph := 0 


( radians ) 


a := 0.0005 


(meters ) 


rlT := 1 10 


(meters ) 


b := a 


(meters ) 


l := 90 


( radians ) 


d := 12.5 a 


(meters ) 


ecd := 1 


( dimensionless ) 


L := 0.5 


(meters ) 


o 

n 

D 


( siemens/meter) 


h = 5 


(meters ) 


nr := 1 5 


(dimensionless ) 
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Define constants: 



o a '= 



0 

w 



(dimensionless ) 



c := 2.9979 10 8 (meters/sec) 



no := 120-31 



1 



C 



f 



(meters/cycle) 




1=0.999 (meters/cycle) 



y := 0 57721 



80 = ■ 



10 * 



9 



36 7i 



(F/m) 



z ' = 0 2 l 



Define angular offset 0: 




2-31 



Define cosine [Ci(x)] and sine [Si ixj ] integrals: 



n := 1 ..50 



Ci(\) -=0.57721 -t-ln(x)- V 



. , . n — 1 2 n 

1 ) \ 

( 2 n) (2 n)' 



Si(\) 



: S 



(- i 

(2 r 



(Q) 



(m 1 ) 



(dimensionless) 



( radians ) 



[ radians ) 



n - 1 i 2 n — I ) 

) N 

- 1 ) (2 n- 1 )' 
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FOLDED DIPOLE ALIGNED ON THE Z-AXIS IN FREE SPACE: 



Equivalent radius a ^: 



ae = exp 



a 2 ln(a) + b 2 ln(b) + 2 a b ln(d) 



(a + b) 



{meters ) 



ae =0.002 



(meters ) 



Electric field E (0) : 



-j k rff 
e J 


cos| 


|~ cos(0)| 


1 - cos l 


!f)| 


J T l u 

2 7i rlT 


sin(0) / 



(V/m) 



En(0) :=— °- 



E|- 



(V/m) 



Minimum Distance to Far-Field r 



mini 



r 0 •= 1 .6 1 (meters ) 



1 1 ;=5L 



(meters ) 



2 (Lr 

r : = (meters) 

K 



irnin *= ma\(r) (meters) 



rrnin = 2.5 



(meters ) 
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Radiation Intensity U ( 0 ) : 



U(0) : = 



rff 

2 T]0 



(|e(0)|) 2 



(W/solid angle) 



Radiated Power P 



radl 



Prad := 2 7i 



U(0)sin(0) dG 



(W) 



Prad = 36.64 



(W) 



Directivity 

Umax := U| — | 



(W/solid angle) 



Do := ^ 71 Umax (dimensionless ) 

Prad 



(dimensionless ) 

Do = 1.641 



Effective Isotropic Radiated Power EIRP: 

EIRP : = Prad Do W ) 

EIRP =60.132 (W) 
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Radiation Resistance R r : 



Rr: = H 
2 * 71 



y+ ln(k-L) - Ci(k l,) + i sin(k L) (Si(2 k L) - 2 Si(k L)> 
2 

+ 1 cos(k L) 4- Ci(2 k L) - 2 Ci(k L) 



(fl) 



Rr = 73.281 



(Q) 



Input Radiation Resistance of linear dipole R ^ : 
Rr 



Rd :=- 



(Q) 



sin 



kl) 



Rd =73.281 



( 0 ) 



Reactance of linear dipole X r : 



Xr:=J12 

4-71 



Si(k L)+ cos(k L) (2 Si(k L,)- Si(2 k L)) 



(Q) 



sin(k-L) 2 Ci(k L) - Ci(2 k D- Ci 



2 k ;ie^ 



Xr =43 142 



<Q) 



Input reactance of linear dipole : 




Xd =43 142 



(Q) 



(Q) 
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Input impedance of linear dipole : 



Zd : = Rd + j Xd 



(Q) 



Zd = 73.281 +43. 142i 



(Q) 



Characteristic impedance of two-wire transmission line Z n : 
Zo = ~ acosh/ — - — ^1 (Q) 

11 UjTb/ 



Zo = 302.312 



( 0 ) 



Characteristic impedance transform for folded dipole Z t : 



Zt:=.j 




(Q) 



Zt =-2.747 10 5 1 



(Q) 



Characteristic impedance of folded dipo Ae_Zmi 

4 Zt Zd 

Zin : = 

2 Zd+Zt 



Zin =293 307 +172.465t 



(Q) 



Input Resistance R s n : 

Rin := Re(Zin) (Q) Rin =293.307 (Q) 
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Voltage Relection Coeffecient I_: 

r ._ Riir-^Zo (dimensionless) V = “0.015 ( dimensionless ) 

Rin -»- Zo 

Reflection Efficiency Cf : 

et .= £cd [ i - ( tr| ) 2 ] (dimensionless ) it = 1 (dimensionless) 



Ga in G i 

G: = itDo ( dimensionless ) GdB = 10 log( G) (dB) 

G = 1.641 (dimensionless) GdB=215 (dB) 



Polarization loss factor PLF: 



=(uir 



PLF 



(dimensionless ) 



Maximum Effective Aperture : 

Aem:=— PLF (m 2 ) 

4 71 



Maximum Effective Height 

hem := 2 | ^ r ^ crn (m) 



VP 



PLF = 1 



(dimensionless ) 



Aem = 0. 1 3 



(m^ 



hem =0.318 



(m) 
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Far Field radiation pattern of vertical folded dipole in free space. 



X(0) := |En(0)| cos 04- 



Y(0) := | En(0)| sin|04-?J 



0 = 0 ° 




0 = 180 ° 
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HORIZONTALLY POSITIONED FOLDED DIPOLE OVER THE EARTH: 



Define angular offsets 01 and 02: 

A i . 71 71 7171 . 

01 = + — — (radians) 

2 2 2 i 2 




Relative complex effective dielectric constant c r p j_ 



a 

crp •= cr- j 

2 71 f co 



component of E-field E (j>: 



E<j»(01 ) : = j rjo - 



-j k rfT 
C J 

no 


cos| 


if 


sin( 01 ) sin(0)J - cosj 


It) 


2 71 1*11 


•J 


^1 - sin(01 )" sin(O)"' 



-J -k-rff 

E4»i =j n°- 

2 ji rfif 




0 component of E-field E 0: 





Z (radians) 



(dimensionless) 



(V/m) 



(V/m) 



(V/m) 



(V/m) 
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(dimensionless ) 



Vertical plane wave reflection coefficient : 



I 2 

I v( 0 j . . _ crp cos(61 )- ^Igp- sin(01 ) 

/ " 2 

op cos(01 ) +■ ^Icrp- sin(01) 



erp cos^02 *) - Jap- sin/02 V 

pvl^ : = - — - — ^ (dimensionless) 

op cos^02_J + Jap - sin (02 J| 2 



Horizontal plane wave reflection coefficient £h • 



1^,(0] ) cos(O) )- Jap- sm( 61 )' 
cos(01 ) -+ Jap - sin(01 



(dimensionless ) 



rhl := 

z 



COS 



COS 



(° 2 ) - Jgp- Sln 

( 02 ) + J°v- sin 




(dimensionless ) 



E-field ( jj> component) with ground plane reflection E ^p : 
E*gp( 01 )- = E«K 0 ])- (e 1 khcos(01, + Th (01 ) e' J khcos(H1) ) (V/m) 




(V/m) 



E 4 »n( e 1 ) •= l E *gP (01 )l 
max( E(j>gp 1 ) 



(V/m) 
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E-field ( 0 component) with ground plane reflection E 0^ p j_ 
E0gp(01) : = E0( 0 1 ) (e 1 k h cos(Gl ’ - I'v(0l ) e‘ J k,vcos(0,) ) (V/m) 



EOgpl, : 



HOI 



( j k h cc 
e 



(te) r , -J k hc<K ( (f2 ) 



- Tvl c 

Z 



(V/m) 



E6n(01) := . l n °PP (01 )l 



max( EOgp 1 ) 



(V/m) 



Radiation Intensity Uh( 0 l): 

Uh(01) :=— [( |E«t>gp(01 )| ) 2 + ( |E0gp(01 )| 
2 no 



(W/solid angle) 



Radiated Power Ph r< ^: 



Phrad -2 n- 



Uh( 0 1 ) sin( 01 ) dOl 



(W) 



Phrad = 55.658 



(W) 



Directivity Dh ^: 



Ugp_max := 



rfi 

2 T|0 



( Eijigp 1 z ) + ( EGgplJ 



(W/solid angle) 



Uhmax := max(Ugp_max) 



(W/solid angle) 
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Directivity Dh n : 



4 -n-Uhmax 
Dho := 

Phrad 



(dimensionless ) 



Dho = 3.704 



(dimensionless ) 



Effective Isotropic Radiated Power EIRPh: 

EIRPh .= Phrad Dho ( W ) 



EIRPh =206.152 



(W) 



Input Resistance Rh ^ n : 



Rhin := 



2 -Phrad 




2 



(Q) 



Rhin = 1 1 1.315 



(Q) 



Voltage Rejection Coeffecient \_ H : 

fH := ( dimensionless ) 

Rhin 4 Zo 



FH =-0.462 (dimensionless) 



Reflection Efficiency th t : 

sht := ccd [ 1 - ( |rii| ) 2 ] (dimensionless ) chi =0.787 (dimensionless ) 
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Gain Gh : 



Gh *= eht Dho 



(dimensionless ) 



GhdB = 10 log(Gh) (dB) 



Gh =2.914 



( dimensionless ) GhdB =4.645 (dB) 



Polarization loss factor PLFv: 



Ehmax:= | E0n( 0ph) | (V/m) 



r 2 h 

Lobes := floorl — + 



oah • = 



0 

0 

Ehmax 



Lobes 



(dimensionless ) 



PLFh := (|o\\ h oah|) 



( dimensionless ) 



PLFh =0.073 



( dimensionless } 



Maximum Effective Aperture Av prn : 

*> 

Ahem ‘= ( ■ ' '■ PLFh (m 2 ) Ahem =0.017 (m 2 ) 

4-h 



Maximum Effective Height hv pm : 

.. . . [Rhin Ahem 

hhem 2- i (m) 



T]0 



hhem =0.142 (m) 
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H-Plane Radiation Pattern (y-z plane) of Horizontal Folded Dipole Over 
The Earth: 



Yl(01 ) := |E<|m(01)| 




X 1 ( 0 1 ) := I H4»n( 0 1 )| COS 





xrcf .XI(Ol) 
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E-Plane Radiation Pattern (x-z plane) of Horizontal Folded Dipole Over 
The Earth: 



Y2(01):= | E0n( 0 1 )| •sin 




X2(01 ) := | E0n( 0 1 )| cos 
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CAGED DIPOLE ANTENNA 
MATHCAD SOFTWARE CAGED_DI.MCD 

This caged dipole application is modeled after a circular array with the 
radiating elements as dipoles. The dipoles are modeled after conductors 
with infinitely small radii and antenna length greater than /./ 4 . Mutual 
interference between adjacent elements are neglected when calculating 
antenna parameters. 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn) . ) 



The caged dipole antenna Mathcad application will compute the following 
parameters : 

VERTICAL CAGED DIPOLE IN FREE SPACE: 

X = Wavelength 

r min = Minimum Distance to the Far-Field 
AF(0,<j>) = Normalized Array Factor 
E ( 0,4> ) = E-Field of Finite Length Dipole 

E t ( 0,<j» ) = Total E-Field 

E tn (M) = Total Normalized E-Field 

U(0.<j>) = Radiation Intensity 

P rac j = Radiated Power 

D Q = Directivity 

EIRP = Effective Isotropic Radiated Power 
R r = Radiation Resistance 

R in = Input Resistance 

I = Voltage Reflection Coefficient 
= Antenna Effeciency 
G = Gam 

PLF = Polarization Loss Factor 
A em = Maximum Effective Aperture 

hem = Maximum Effective Height 

flow = Low Frequency 
fhigh = High Frequency 
BW = Bandwidth 

VERTICAL CAGED DIPOLE OVER EARTH : 

c rp = Relative Complex Effective Dielectric Constant 
AFv(01,<t)) = Normalized Array Factor 
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E v ( 01 , 4> ) = E-Field of Finite Length Dipole 

T v (Gl) = Vertical Plane Wave Reflection Coefficient 

Ev t (01 , <(i) = Total E- Field 

Ev tn {0l,$) = Total Normalized E-Field 

Uv ( 01 , 4» ) = Radiation Intensity 

Pv ra d = Radiated Power 
Dv 0 = Directivity 

EIRPv = Effective Isotropic Radiated Power 
Rv r = Radiation Resistance 
Rv in = Input Resistance 

rv = Voltage Reflection Coefficient 
cv t = Antenna Efficiency 
Gv = Gain 

PLFv = Polarization Loss Factor 
Av em = Maximum Effective Aperture 

hv e m = Maximum Effective Height 

HORIZONTAL CAGED DIPOLE OVER EARTH: 

AFh(01 / (f)l) = Normalized Array Factor 

(01 ) = Horizontal Plane Wave Reflection Coefficient 

E0h(01) = 0 Component of E-Field 

E0h t (01 , <J)1 ) = 0 Component of Total E-Field 

E0h tn (01 , <(>1 ) = 0 Component of Total Normalized E-Field 

E((>h(01) = <(> Component of E-Field 

E(j)h t (01 , (j)l } = Component of Total E-Field 

E^hfn ) = $ Com P onent of Total Normalized E-Field 

Uh (01 r (jil ) = Radiation Intensity 

Ph ra d = Radiated Power 
Dh Q = Directivity 

EIRPh = Effective Isotropic Radiated Power 
Rh r = Radiation Resistance 

Rh^ n = Input Resistance 

TH = Voltage Reflection Coefficient 
eh t = Antenna Efficiency 
Gh = Gain 

PLFh = Polarization Loss Factor 
Ah e m = Maximum Effective Aperture 
hh em = Maximum Effective Height 
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The following known or estimated data must be entered: 



N 

a 

L 

h 

f 



W w 

0Vp 

ah 

0 h 



w 



P 

j = 



Number of Dipoles 
Radius of Antenna 
Length of Antenna 
Antenna Height Above Ground 
Frequency of Interest 
= Antenna Feed Current 

= Input Impedance 

= Relative Dielectric Constant of Ground 
Conductivity of Ground 

= Incoming Wave E-Field Unit Vector for Antenna (Free Space) 

= Antenna Unit Polarization Vector (Free Space) 

= Incoming Wave E-Field Unit Vector for Vertical Antenna 

= Direction of Incoming Unit Vector (Vertical Antenna) 

= Incoming Wave E-Field Unit Vector for Horizontal Antenna 

= Direction of Incoming Unit Vector (Horizontal Antenna) 

Number of Increments in Degrees for Far Field Radiation Pattern 



Enter input data here: 



GW • = 


fl 

0 


(dimensionless ) 


aa ’= 


1 

0 




\o) 






A 



cw := 



ohw 



'o' 

0 



( dimensionless ) 



Ovp 



80 n 
180 



(dimensionless ) 



0hp : = 0 



(dimensionless ) 



( dimensionless ) 



(dimensionless ) 



N := 8 



( elements ) 



Io := 1 0 (A) 



1 

a -= - (meters ) 

8 



Zo • = 75 0 (Q) 
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L := 0 5 



(meters ) 



j := 90 



( increments ) 



h .= 0.25 



(meters ) 



o .= 0.01 



f := 3 10 8 (Hz) 



a* := 1 5 



Define constants and calculate wavelength: 
c := 2.9979 10 8 (meters/sec) rjo := 1 20 71 



l := 



c 

f 



(meter s/cycle ) 



k := 




/. 




co 




36ti 



(A) 



( farads/m) 



rtr •= 10 5 



ccd := 1 



l =0.999 



(meters/cycle ) 



z -=0 2 j 



Calculate caged dipole parameters parameters: 
Define angular offset 0_: 

0 " = 0 . — — - . 2 * 7i (radians) 4> ' = 0 , — 2 n 

j j 



( siemens/m) 



(dimensionless ) 



(Q) 



(m !) 



(meters ) 



(dimensionless ) 



( increments ) 



( radians ) 
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Minimum Distance to Far-Field r m . n : 
r Q := 1 .6 /. (meters ) 



r, : =5i/i 



L 2 4 ( 2-a) 2 



(meters ) 



2 [l 2 + (2 a) 2 ] 
X 



(meters ) 



rmin — max(r) (meters) 



nmn = 2 795 



VERTICAL CAGED DIPOLE IN FREE SPACE: 



Array Factor AF( 0 ): 



n:= 1 N 



( elements ) 



d) := 2 n — 

n N 



AF(0,<j>) := 'y ^ In e 



j kasin(0) •cos (♦- '0 



E- Field of finite dipole E ( 0 ) : 



E(0) :=j r|o 



In e 



- j k rfi 



/k L . . 

CO.S I — — cos(0)| - cos 



(W 



71 ill \ 



sm( 0) 



(meters ) 



( radians ) 



(dimensionless ) 



(V/m) 
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(V/m) 



Total E-Field E ,-( Q)& Normalized E-Field E i- n (6): 

Et(G, <(>) •= E(0) 'AF(M) 



Etn( 0, <J>) 



Et( 0 , 4>) 




(V/m) 



Radiation Intensity U ( 0 ) : 

U(0,f) :=— (|Et(0J)|) 2 
2-r|0 



(W/solid angle) 



Radiated Power P rp ,H : 



2 n 



Prad •= 



71 

U(0,<j>)sin(O) d0 d<t> 



Jo Jo 



(W) 



Prad = 0.44 1 



(W) 



Directivity D n : 



Uma\:=U| — ,0[ (W/solid angle) 



Do = 1.543 



Do := 



4 7i Umax 
Prad 



(dimensionless ) 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 

EIRP := Prad Do (W) EIRP =0.681 (W) 
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Radiation Resistance R r : 

2 Prad 

Rr : = -- (Q) Rr = 56.487 (Q) 

(|ln |) 2 



Input Resistance R ln : 

Rin •= — (fi) 




Voltage Rejection 

r;= Rin-Zo (d 

Rin + Zo 



Coeffecient Y: 



imensionless ) 



Rin = 56.487 (Q) 



I =-0141 (dimensionless ) 



Reflection Efficiency : 

ct •= ccd [ 1 - ( |r| ) 2 ] (dimensionless ) a = 0.9S ( dimensionless ) 



Ga in G i 

G-^etDo ( dimensionless ) GdB •= 10 Iog( G) (dB) 

G = 1.513 (dimensionless ) GdB = 1 79S (dB) 



Polarization Loss Factor PLF : 



4U»I) 2 



PLF 



(dimensionless ) 



PLF = 1 



(dimensionless) 
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Maximum Effective Aperture : 



C' \ 

Aem := 1 ' PLF (m^) Acm =012 (m^) 

4 71 



Maximum Effective Height 

hem =2- j ^ cm (m) hem =0.268 (m) 

N Tio 

Define Frequency Relationship F(x): 

F(x) 5.696 10' 5 \ + 0.3 1481 (dimensionless ) 



Frequency Range f low and fhigh : 






flow := f- p|— j f (Hz) 


flow =2.056 10 8 


(HZ) 


fhigh =f4F|— ] f (Hz) 


thigh =3 944 10 8 


(Hz) 


\ 2 a/ 







Bandwidth 


BW: 








BW := thigh - 


flow 


(HZ) 


BW = 1.888- 10 8 


(HZ 



154 



E-FIELD PATTERN. This far-field radiation pattern is for the vertical 
caged dipole in free space: 



Number of Dipole Conductors: 



Radius of Antenna: 



X(0) := |Etn(0,<|>)| •cost 0 — : 



N =8 

a =0.125 



(meters ) 



Y(0) := |Etn(0.(jp)| sinlO- - 



0 = 0 ° 




0 = 180 ° 
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VERTICAL CAGED DIPOLE OVER EARTH: 



Define angular offset 0_ 1 and j) 2 : 



A! 71 71 71 71 , , . 

01 (radians) 

2 2 2 j 2 



<t>2 - — z 

Z J 



02 = --4— z 

z 2 2j 



( radians ) 



Relative complex effective dielectric constant c r p : 



erp = er- j 



2-Tifeo 

Array Factor AFv( O ): 

V I k a sin< 01 ) cos 

Ine v 0 

' ' Z j 

n 

V“I T J k a s,n(02^ cos^-Oq) 

AFvgp z = > III e 
n 

E- Field of finite dipole Ev( Q l): 



Ev(01) =j r| 0 - 



t -j kriT 

In e * 


^cos| 


If cos(GI )j 


| - cos| 




2 H 1 11' 


\ sin(01) J 



Evgp: = j k|0 



In e 



-j k iff 



2 • 71 I'll' 



cosf— cos/02 Vl - cosl—'j 

1 2 * v \ 2} 



i* 2 ) 



Vertical plane wave reflection coefficient f v (Ql ) : 



rv(01 ) : = 



ap cos(01 ) - J 


Lp- sin(01 y 


op cos(01 )*k 


U- sin (^i y 



( radians ) 



( dimensionless ) 



( dimensionless ) 



( dimensionless ) 



(V/m) 



(V/m) 



(dimensionless ) 
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(dimensionless ) 



erp cos| 


N - J 


ferp- sin | 




op cos| 


M 




fop- sin | 


M 2 



Total E-Field EV t-(G l) & Normalized E-Field Ev t - n (6 l): 

Evt(eM) = Ev( 01 )-AFv (01 ,4.) (e 1 k h cos< 01 ' 4 I v( 01 ) e' J K h cos( 01 ’) (V/m) 



E\lgp z := 



c AC { i khcos ( te ) r 1 *J khcos ( e2 )\ 

Evgp z AFvgp^ I e + fvl e I 



(V/m) 



Evtn( 01 , 4 >) 



|E\t(ei,»)| 

max(Evtgp) 



(V/m) 



Radiation Intensity Uv( 6[) 
Uv(01, <(,):= — (| Evt(6M)|) 2 

2 T]0 



(W/solid angle 



Uv 






2 rjo 



; W /solid angle) 



Radiated Power Pv r ^rj : 



P\Tad •- 



Uv(01 . <j>) sin(01 ) dO 1 d(J> 



(W, 



Pvrad = 0.168 



(W) 



Directivity Dv n : 

Uvmax Uv(0vp,O) (W/solid angle) 
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Directivity Dv n : 



~ 4 n Uvmax 

Dvo := 

Pvrad 



Dvo = 2.803 



Effective Isotropic Radiated Power EIRPv: 

EIRPv := Pvrad Dvo (W) EIRPv = 0.47 1 



Radiation Resistance Rv 



r_i 



Rvr := 



2 Pvrad 
( I In I ) 2 



IQ) 



Input Resistance Rv ^ n : 
Rvr 



Rvin *=- 



(Q) 



. (k L 

sin 

2 



Rvr =21.515 



Rvin =21.515 



Voltage Relection Coeffecient I V: 

rV := ^ Un — — (dimensionless) I'V =-0.554 

Rvin +- Zo 



Reflection Efficiency zv t : 

evt : = ccd [ 1 - ( |rv| )"] (dimensionless) cvt =0.693 



(dimensionless ) 



( dimensionless ) 



(W) 



(Q) 



(Q) 



(dimensionless ) 



(dimensionless ) 
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Gain Gv: 



Gv := evt Dvo 



(dimensionless ) 



GvdB - lOlog(Gv) (dB) 



Gv = 1.942 



(dimensionless) GvdB = 2.883 (dB) 



Polarization Loss Factor PLFv: 
Evmax = Evtn( 0 vp,O) (V/m) 



Evmax 



ova * = 



lobesv 

0 

0 



(dimensionless ) 



lobesv = floor 2 - + 1 



PLFv := ( |owova|) 



(dimensionless ) 



Maximum Effective Aperture 
2 

A vem := i£L_L PLF v ( m 2 ) 

4 71 



Maximum Effective Height hv prn : 



, . [RvrAvem 

hvem .= 2 / — — (m) 

1 T]0 



PLFv = 0.689 



(dimensionless > 



Avem =0 106 



(m 2 ) 



hvem =0 156 



im) 
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E-FIELD PATTERN: This far-field radiation pattern is for the vertical 

caged dipole over earth: 



Number 


of 


Dipole Conductors: 


N =8 


( elements ) 


Radius 


of 


Antenna : 


a =0.125 


(meters ) 


Height 


of 


Antenna Above Earth: 


h = 0.25 


(meters ) 



Xl(01) := |Evtn(ei,*)| coslOl + — 



Y 1(01 ) := [Evtn(01 ,<{.)[ sin 01 +- 



0 = 0 ° 




xief m .Xl(ei) 
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HORIZONTAL CAGED DIPOLE OVER EARTH: 

Define angular offset <ji l and _j>2 : 

4>1 -= 0 2* 7t (radians) 

J 



Array Factor AFh ( 9 1, $ 1 ) : 



n .= I N 



(elements ) 



0 .= 2nl 
n N 



AFh(91,*l) =2 In e 



ka (sin( 01 ) sin( 01 ) sinfir'^ cos( 01 ) -cosf©^ 



AFhl ? :=^In e 
n 



j k a (sin^O^I sin^02^l sm^6^j -+cosl^02^l cos^O^ljl 



0 component of E-Field of finite dipole E O h( O l): 



t -i krfl 

E0h( 01) •= j T|o — il 

2 /i rff 



cos|^— sin( G 1 ) | - cos 



IT- 



y 1 - sin(01 )" 




T -j k rlT 

EOhl •= j rp — 

2 2 rtf 




Total E-Field E O h ( 0 1, & 1 ) : 

E0ht( 01 , 4 .] ) := E0h(01 ) AFh(01 ,4-1 ) (e' khcos(ei) - Iv( 01 ) e‘ J k h cos(H1) ) 



E0htl = 

Z 



E0hl • AFhl 



/ 1 k h cosfo'l -] 

w - rvl ? e 



k h e os (02^1 ) 



( radians ) 



(dimensionless ) 



(dimensionless ) 



(V/m) 



(V/m) 



(V/m) 



(V/m) 
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(V/m) 



Normalized E-Field E Oh t n (0 1 , jt l ) : 

E0hm(0l,*l) = J2W0MI)| 
max(E0htl ) 



Horizontal plane wave reflection coefficient £>^(01) 



rh(01) : = 



cos(01 ) - aJ 


^irp - sin(01 ) 2 


cos(01 ) 4 J 


Irp- sin(01 ) 2 



Hil = 



cos| 


w 


'-■J 


[erp — sin| 




cos| 




[erp- sin| 





component of E-Field of finite dipole E<j>h(0l) 



H<t>h(01 ) — j -r]o 



In e 



-j k rfT 



2 7i rff 



1 - cos 



i 



C1U1 . In e' j k rff ( ( k L^\ 

2 ' 2-wrff \ l 2 // 



Total E-Field E $h t (fl l , & 1 ) & Normalized E-Field E ^h,- n (0: 

E4»ht(01 ,<M ) = E4*h(01 ) AFh(01 ,<M ) (e' k h cos(6l) + ni(61 ) e' J khcos<ei) ) 



E^htl := 



/ j khcos(()2J -) k h cos(02j \ 

E<t>hl z AFhl^e k * -f Hi 1^,0 w l 



E<t>htn( 01 , 4»1 ) zjEMeMOl 
max(E<i>htl ) 



(dimensionless) 



(dimensionless) 



(V/m) 



(V/m) 



, 4> 1 ) : 



(V/m) 



(V/m) 



(V/m) 
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Radiation Intensity Uh( Q l, $ l): 



Uh(01 ,$1 ) •=— [( |E0ht(Ol ,4.1 )| ) 2 4 ( 1 E<J»ht(0 1 .$1 )| ) 2 ] 

2r]0 



(W/solid angle) 




(W/solid angle) 



Radiated Power Pv r ^ H : 



2-71 



71 

0 



Phrad •= 



Jo 



Uh(01 , <() 1 ) sm(01 ) dO 1 1 

<0 



(W) 



Phrad = 0.739 



(W) 



Directivity Dh n : 

Uhmax •= Uh(0hp,O) (W/solid angle) 



Dho 



7 i' Uhmax 



Phrad 



( dimensionless ) 



Dho = 4.662 



(dimensionless ) 



Effective Isotropic Radiated Power EIRPh: 

EIRPh := Phrad Dho (W) EIRPh = 3 443 (W) 



Radiation Resistance Rv r : 

0 . PhrnH 

Rhr •= — (Q) Rhr = 9-1532 (Q) 

( I In | ) 2 
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Input Resistance Rh -j n : 



Rhin := — — — (Q) 




Voltage Relection Coef fecient flL 

rH := ^ in — — (dimensionless) 

Rhin-t Zo 



Rhin =94.533 (Q) 



FH =0.115 (dimensionless) 



Reflection Efficiency gVf : 

eht * = ecd [ 1 - ( | rH| ) 2 ] (dimens ionless) chi =0.987 (dimensionless) 



Gain Gh i 

Gh cht Dho (dimensionless) *• GhdB '= 10 log(Gh) (dB) 

Gh=4.6 (dimensionless) GhdB =6.627 (dB) 



Polarization Loss Factor PLFh : 



Ehmax * = E<|)htn( 0hp , 0 ) ( V / m ) 



aha • = 



0 

0 

Ehmax 



lobesh 



( dimensionless ) 



PLFh 



• = ( |ah\v cha|) 



(dimensionless ) 



lobesh '= floor [2 — + 1 



PLFh =0.958 ( dimensionless ) 
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Maximum Effective Aperture Ah ^ : 



Ch > 2 

Ahem := — — - PLFh (m 2 ) 

4-a 



Ahem =0.35 



(m 2 ) 



Maximum Effective Height hh ^ : 



,, _ Rhr Ahem , 

hhem =2 J - (m) 

J T]0 



hhem = 0.593 



(m) 
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H-PLANE RADIATION PATTERN: This far- field radiation pattern in the 

(y-z) plane is for the horizontal caged dipole over earth: 



Number 


of 


Dipole Conductors: 


N =8 


(elements ) 


Radius 


of 


Antenna : 


a =01 25 


(meters) 


Height 


of 


Antenna Above Earth: 


h = 0.25 


(meters ) 



X2(01 ) := | E4»htn( 0 1 ,0)| cos 




Y2(01 ) := |E<t>htn(01 ,0)| sinfoi + 



0 = 0 ° 
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E-PLANE RADIATION PATTERN: This far-field radiation pattern in the 

(x-z) plane is for the horizontal caged dipole over earth: 



Number of Dipole Conductors: 
Radius of Antenna: 

Height of Antenna Above Earth: 



N = 8 (elements ) 

a = 0 125 (meters ) 

h = 0.25 (meters) 
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THE HELICAL ANTENNA (REFLECTOR OPTION) 

MATHCAD SOFTWARE- DISH^HEL . MCD 

When built to the proper specifications, the helical antenna possesses 
many qualities which make it suitable for a wide variety of 
communications applications. If the following conditions are satisfied 
the helix will exhibit a highly directional axial main lobe, low side 
lobe level, negligible mutual interference with adjacent antennas, low 
voltage standing wave ratio (VSWR) , and resistive input impedance over a 
wide frequency band: 

• 8 1 C A 1 1 * 15 

n > 3 

12 < a < 14 

(Note: X in a subscript indicates the dimension is in wavelengths. 

Mathcad equations can not use symbolic subscripts. Therefore, the symbol 
X will immediately follow the parameter m equations (i.e., Ck) m lieu 
of subscripts.) 

The helical antenna Mathcad application will compute the following 
parameters (Items with * indicate parameters that are calculated for 
both axial and peripheral feed geometries): 

C = Circumference of Helix 
X = Wavelength 
a = Pitch Angle 
Do = Directivity 
p = Relative Phase Velocity 
Eq ^ = Electric Field Components 

4 ; = Array Factor Phase Shift 

U = Radiation Intensity 
Prad = Radiated Power 
R = Antenna Input Resistance* 

T = Voltage Reflection Coefficient* 

e r = Reflection Efficiency* 

h em = Maximum Effective Height* 

G = Gain* 

EIRP = Effective Isotropic Radiated Power* 

A em = Maximum Effective Aperture* 

AR = Axial Ratio 

PLF = Polarization Loss Factor 

BW = Bandwidth 

fhigh = Upper Frequency Limit 
flow = Lower Frequency Limit 
Acceptable Conductor Diameter 

E x,y,z = Electric Field Cartesian Coordinates 
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0p/(j)p = Unit Polarization Vector Coordinate Angles 

o a = Antenna Unit Polarization Vector 
r m in = Mimnum Distance to the Far-Field 



The following data must be input based on known or estimated data: 

D = Diameter of Helix (Center to Center) 

S = Spacing Between Turns (Center to Center) 

L = Length Along Conductor of One Turn 
n = Number of Turns 
d = Diameter of Helical Conductor 
f = Frequency of Interest 
m = Desired Mode 

1 = Number of Increments in Degrees for Far Field Radiation Pattern 
I Q = Antenna Feed Current 

Z Q = Cha racten Stic Feed Impedance 

o w = Incoming Wave Electric Field Unit Vector 



Enter input data here: 




aw = 



(dimensionless . 



0 



D := 



0.1 



L := 0.097 (meters) 



(meters ) 



n 



n := 10 



( turns ) 



f := 3 10 10 (Hz) 



(dimensionless ) 




(Q) 
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S := 0.023 
d =0.005 

x := 1 

y 



(meters ) 
(meters ) 

(meters ) 

(meters ) 



l = 360 (degrees ) 

lo = 1 (A) 

z ■= 1000 (meters ) 



Calculate helical geometric parameters and define constants: 

.8 



c •= 2.9979 10 (meters/sec) 



r| 0 -= 1 20 Ti (Q) 



C ‘= n-D (meters ) 






(meters/cycle ) 



C = 0. 1 (meters ) 



).= 9.993" 10“ (meters /cycle ) 



C). '•= — 



(dimensionless ) 



a := atan — (radians) 

\C 



C/. = 10.007 (dimensionless ) 



a =0 22607 (radians) 



T* - L 

La • = - 



S 

(dimensionless) S/. • = - (dimensionless) 



La =9 70679 (dimensionless) SI = 2.30161 (dimensionless) 



ad •= -a (degrees) 



k := 2 - 



(m 



ad = 12 95276 (degrees) 
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Define angular offset 8 
from helical axis: 



Minimum Distance to the Far-Field 



£minl 



i:=0,-..2n 



( radians ) 




( radians ) 



= 1 . 6 /. 



Tj := 5 n S 




rmin '= max( r) 



rmin = 10.5874 1 



Relative Phase Velocity c: 



L/. 

p := 

S/. 4- m -4* 




(meters ) 



(meters ) 



(meters ) 



(meters ) 



(meters ) 



{ dimensionless ) 



p =2 89616 



(dimensionless) 



Electric Field Components E ft , E ^: 



E(0) ■- 



sm (rJ 


sinj 


f n i)i(0)j 


sinj 


ifi: 


1 



cos(G) 



(V/m) 



E0(0) :=E(0) 



(V/m) 



Em := j E(0) 



(V/m) 



Radiation Intensity U ( 0) 

.2 



U(0) •=— •( |E(0)| ) 2 

% 



(W/solid angle) 



Radiated Power P 



rad-L 



Prad := 



Z • 71 



U(0) sin(0)dG cty 



(W) 



Prad = 1 67351- 10” 3 



(W) 



Directivity D n : 



2 

Do ‘=12 Cl n SI (dimensionless) 



Do2 '= - — (dimensionless ) 
Prad 



Do = 2 7658 10 (dimensionless ) Do2 = 19.91822 (dimensionless ) 
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Axial Ratio AR : 



AR := 



Lb |sin(a) - — 



AR = 1.17586 

Effective Isotropic Radiated Power EIRP: 

EIRP;=Prad Do (W) E1RP2 : = Prad Do2 



EIRP =46.286 (W) 



EIRP2 =0.03333 



Polarization Loss Factor PLF : 



Bp - a tan 



r~i \ 

■yx -i-y 



( radians ) 



4>p = atan - 

\x 



Bp = 1 41421 10 (radians) 



d>p =0 7854 



Ex =EB(Bp) cos(0p) cos(<t>p)- E<KBp) sin($p) 



Ey =EB(Bp) cos(Bp)-sin(<)>p) + E<KBp) cos(^p) 



Ez =EB(Bp) sin(Bp) - 1 



oa =- 



J(|Ex| ) 2 +- ( |Ey| ) 2 + ( |Ez| ) 2 



lx' 

iEz f 



oa = 



0.5- 0. 5j 
0.5 + 0.5j 

-9 99999M0 



“4 



(dimensionless ) 



(dimensionless ) 



(W) 



(W) 



( radians ) 



( radians ) 



(V/m) 



(V/m) 



(V/m) 



(dimensionless ) 



(dimensionless ) 
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Polarization Loss Factor PLF: 




(dimensionless ) 



PLF = 1 



Radiation Resistance R r : 

Rr := 2 Prad (Q) Rr = 3.34702 10~ 3 

(IH ) 2 



Dual Parameters 



Axial Feed 
Input Resistance R: 
Ra:= 140.^0. (Q) 

Ra = 4 42874 10 2 (Q) 



Peripheral Feed 



Rp := 



150 

Ja 



Rp =47 41756 



Voltage Reflection Coef f ecient J : 

Ra- 



Ta . 



o 



Ran- Z, 



Ta =0 49399 



(dimensionless ) 



(dimensionless ) 



Rp-Z 0 
Ip ■= 

R P- Z o 

Ip =-0.51962 



Reflection Efficiency c r : 

era .’= 1 - ( | fa | ) 2 (dimensionless ) 



txp := 1 - ( |rp| ) 2 



£ra = 0.75597 



(dimensionless ) 



crp = 0 72999 



(dimensionless ) 



(Q) 



(Q) 



(Q) 



(dimensionless ) 



(dimensionless ) 



(dimensionless) 



(dimensionless ) 
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Ga in G i 



Ga := era Do (dimensionless) Gp := op Do (dimensionless ) 



Gadb:= lOIog(eraDo) (dB) 



Gpdb := 10 log( ap Do) (dB) 



Ga = 2.09088' 10 4 (dimensionless) Gp = 2.01902* 10 4 (dimensionless ) 



Gadb =43 20328 (dB) 



Gpdb =43.05 14 (dB) 



Ga2 •= era Do2 (dimensionless) Gp2- = erpDo2 (dimensionless ) 



Gadb2 := 10 log(era Do2) (dB) 



GpdbZ ■= 10 log(aj> Do2) (dB) 



Ga2 = 15 05765 ( dimensionless ) Gp2 = 14 54017 (dimensionless ) 



Gadb2 = 1 1 77757 (dB) 



Gpdb2 = 1 1.62569 (dB) 



Maximum Effective Aperture A, 



emu 



A . era ). Do r 
Aema •= PLF 

4 n 



(m 2 ) 



A ap /. Do nT y 

Aemp = — PLF (rrr; 

4 n 



Aema =0.16615 



(m 2 } 



Aemp =0 16044 



(m 2 ) 



A ~ . era >. 2 Do2 nT 
Aema2.= PLF 

4 71 



(m 2 ) 



. 2 . 



. ^ . txp /. Do2 o . 

Aemp2 •= — PLF (m z ) 

4 7t 



Aema2 = 1 19657 10 



(m 2 ) 



Aemp2 = 1.15545 10 4 (m 2 ) 
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Maximum Effective Height h ^ m : 
hema :r 2 Irt (m) 

J % 

hema = 2.42912" 10 3 (m) 

hema2 := 2 |Rr ^ em . a ^ (m) 

J % 

hema2 =6.51873 10” 5 2 (m) 



hemp := 2- Rr 



Aemp 



hemp = 2.38701" 10 



r 3 



/pr. Aem P 2 



hemp2:=2 Rr 



hemp2 =6.40574 10 



rs 



Bandwidth : 

, 1.15C 

thigh:- (Hz) 

C 

thigh = 3.44759 10 9 (Hz) 

BW := thigh - flow (Hz) 

BW = 1.04927 10 9 (Hz) 



flow 



.8 c 

~ 



flow =2.39832- 10 9 



Acceptable Conductor Diameter: 

dmin := .005 >. (m) dmax:= 057. 



dmin = 4 9965 10 5 



(m) 



dmax = 4 9965 10 



(m) 



(m) 



(m) 



(m) 



(Hz) 



(Hz) 



(m) 



(m) 
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HELICAL ANTENNA FAR-FIELD RADIATION PATTERN 



For the purpose of this far-field radiation pattern, the helical antenna 
axis is equivalent to the Ex = 0 grid line. The pattern is essentially 
symetric when rotated about the antenna’s axis. 



Ex(0) := E 




cos(0) 



Ey( 6) *=E|0-- 



•sin(0) 




— E FIELD 



Ex< 0 ) 
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PARABOLIC REFLECTOR ANALYSIS WITH HELICAL ANTENNA FEED 



The helical antenna Mathcad application (Reflector Option) will compute 
the following additional parameters (Note: a 2 or 22 will follow most 
secondary field parameters in this application, and inorder to speed up 
computation time corresponding vectors are created with subscript A & B 

P x ^y = E-Field Integral and Summation Equation Vector Components 
E aX/ ay = Aperture E-Field Components 

0*,p = E-Field Integral and Summation Equation Parameters 

Eq = Secondary Field Theta E-Field Component 

E^ = Secondary Field Phi E-Field Component 

U = Secondary Field Radiation Intensity 

?rad = Secondary Field Radiated Power 

D q = Secondary Field Directivity 

G = Secondary Field Gain 

8 ap = Aperture Efficiency 

c s = Spillover Efficiency 
= Taper Efficiency 
Cp = Phase Efficiency 
c x = Polarization Efficiency 
c ohmic = Ohmic Efficiency 
A e = Reflector Effective Area 

h e = Reflector Effective Height 

EIRP = Reflector Effective Radiated Power 
E X/ y fZ = Secondary Field Cartesian Components 

o a = Reflector Unit Polarization Vector 

PLF = Reflector Polarization Loss Factor 

r min = Minimum Distance to Secondary Far-Field 

0 O = Reflector Periphery Coaltitude 

0p,<t>p = Theta and Phi for Polarization Loss 



The following additional data must be inputted based on known or 
estimated data : 

E Q = Electric Field Scale Factor 

a ref = Radius of the Mouth of the Reflector 

foe = Reflector Focal Length 

£b = Blockage Efficiency 

e S p = Spar Efficiency 

tl,t2 = Radiated Power Increments for Phi and Theta 
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il = Secondary Field Increments 
N, M = Summation Increments 

r f f 2 = Secondary Field Observation Distance 

x,y,z = Polarization Loss Factor Coordinates 
ct w = Incoming Wave Unit Polarization Vector 



Enter additional input data here: 

foe := 5 (meters) aref:=5 

(Note: (f) should be greater 

than minimum distance to the 
far-field for the feed antenna) 



cb r .96 



esp = 95 



(dimensionless ) tl =10 



12 = 10 



( increments ) 



il = 10 



N = 10 



(increments ) 



M = 10 



rfn = 10 ? 



(meters) 



x = 1000 



cw 2 



0 



> =1000 



(dimensionless ) 



z = 1000 



(meters ) 



( dimensionless ) 



( increments ) 



( increments ) 



( increments ) 



(m) 



(m) 



(m) 
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Establish Integral and Summation Increments: 






<H :=0,— ..2-ti 
tl 



( radians ) 



A := 1 . N 



*2 := 0,— .. 2 it 

il 



( radians ) 



B := 1 M 



02 






icr 6 , — i a 6 ..-- io ' 6 

1118 36 36 



rp := 0,-^f . aref (m) 

t2 



^ :=2 n ~ 



M 



rp 1 := aref— 

A N 



(m) 



922 := I- --- 
A \ 1 8/ N 36 



02 2„ : = -• 



36 



Calculate reflector geometric parameters: 



0o := 2-atanf-"— ^ 

U foe } 


( radians ) 


r o : =16 7 


00 = 0.9273 


( radians ) 


i j := 10 aref 


4 foc 2 + rp" 
4 roc 


(m) 


8 aref 2 
^ JL 


0p(rp) := 2 alanf-— | 
\2 foe/ 


( radians ) 


mun := max( r) 


(N 

- 

J* 

^r 

ii 

< 

c. 


(m) 


rmin = 2.0014 




( radians ) 





( increments ) 



( increments ) 



( radians ) 



( radians ) 



( radians ) 



( radians ) 



(m) 



(m) 



(m) 



(m) 



(m) 
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Aperture Electric Fields E a XJ , E a y j_ 

Eax(rp,4>l ) := (E0(0p(rp))--cos(^l ) + E4>( 0p( rp) ) sin( 4>1 )) e* k ( p(rp) “ foc 2) (V/m) 

Eav(rp,4>l ) = (E0(0p(rp)) sm( 1 )4- E<H9p(rp)) cos(<t>l )) e 1 ' k(p(rp) ~ fc>w) (V/m) 

Eaxl A,B (-C 05 ^ 1 ! b )) + E^OplJ si n (<bl ] e k ( rl A- foc -) (V/m) 

Eayl A B := (Ee(0p]^ sin^l + E^Opl^ cos^l 1 B )) e k ^' A ^ (V/m) 



Electric Field Integral Equation Vector Components P v , v j_ 
"aref 



Ax(02,*2,(|»l ) := 



Ea\(rp 4>1 ) c' k (cosiel) cos(«2)-fsin(iH) sm < ll,2 )). r p C (jp 



0 



(V/m) 



Px( 02 , <j>2 ) := 



Ax(G2 .4.2,^] ) d(M 



(V/m) 



Ay(92,4>2,4>] ) := 



aref 



Eay(rp,*])e' ^"(^Mcosiei) cos<*2 ) + s,n<?i).s,n(«2)) ^ ^ (V/m) 



P> ( 92 , 4»2 ) '• = 



Av<02,*2.<M ) d4>l 



(V/m) 



Axl 



A,B 



aref 

0 

*2n 



Eaxl e J (cos(oil a \ cos(e22^^sm(cll^ sinM) (v/m) 

A, B 



Pxl 



A.B 



A - 1 A.B ' = 



Axl AB dpll 



(V/m) 



0 



Pyl 



A,B 



aref 

0 

2 n 



Eay ] C J k cos ( e22 ^ +S "’( CI 1 B I “W) rp, drp, (V/m) 



A >' ] A.B d611 



(V/m) 



0 
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Secondary Field Electric Field Components E g, 

• | -j k rfT2 

. J -k e J 



E02(02 , 4*2 ) •=- 



2 7i rfl2 



(Px( 92,^2) cos(^2)+ P>( 02,4,2) sin($2)) 



i -j k riT2 

E<j)2(02 ,d>2) .= i— ( Px( 02 , 4»2 ) - sin(<t>2) + P}^(02^2) cos(<t)2))-cos(02) 

2-71 rff2 



• i, -j krfT2 

e93 a,b 2 J (^'a.b “*(* 22 b) * p > 'a.b s '"(* 22 b)) 



i .j *k-rfT2 

e * 3 a,b = J 2 Irfn [ Px1 a.b (- sin (* 22 B )) 4 cos (* 22 d ] cos ( e22 J 



Secondary Field Radiation Intensity U: 

U2( 02,^2) := — ■[( | E02 ( 02 , <J»2 ) ( ) 2 + ( |E*2(02 ,*2)| ) 2 ] 

2 % 



U3 a.b™[(I e 83 a.bI)M| E * 3 a.bI) : ] 

^0 



Secondary Field Radiated Power P rPi H : 

Prad2 (yrl (i^) -ZE^K’* 2 ^ 



Prad2 =0.19318 



Secondary Field Directivity : 

Dq22 4 nmax(U3) 

Prad2 

Do22 = 1.04762 10 5 



(v/m) 



(v/m) 



(V/m) 



(v/m) 



(W/solid ang) 



(W/solid ang) 



(W) 



(W) 



(dimensionless ) 



(dimensionless ) 
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Spillover Efficiency : 
’00 

U(0p) srn(Op) dOp 

•0 

ES ' = 

*7t 

U(0p) sin(Gp) dOp 

• 0 



es =0.19013 



Taper Efficiency c t : 
f Go 



Et '= 



Pgr “i 1 * 



0 



2 cot 



rGo 



U(Gp) 4 n 
Prad 



sinf Op ) dOp 




et =0 50837 



Phase Efficiency : 

-4x6 10 ~ aref 

Ep •= e 



Ep = 0.68574 



Aperture Efficiency c a p j_ 

ex - 98 (dimensionless) cohnnc •- 98 



cap es d Ep ex cb csp cohmic (dimensionless cap =0 05805 
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(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP22 := Prad2Do22 



(W) 



EIRP22 = 2 023 7 9 1 0 4 (W) 

Ga in G i 

G22 — cap Do22 (dimensionless ) G22dB := 101og( G22 ) (dimensionless ) 

G22 =6.08192 10 3 (dB) G22dB = 37.84041 (dB) 



Effective Aperture : 



Ae22 cap ti aref 



(m 2 ) 



Ae22 =4.55961 



(m 2 ) 



Radiation Resistance R r : 



Rr22 



Prad2 

(|lo |) 2 



( 0 ) 



Rr22 =0.38636 



(Q) 



Effective Height : 



J Ae22 

Rj 22 — — (m) 

% 

he22 =0.13672 (m) 
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Polarization Loss Factor PLF : 



ep: 




( radians ) 



*p2 




0p2 =0.95532 



( radians ) 



4>p2 =0.7854 



E04 := E02( 0p2 , <t>p2 ) (V/m) 



E<{4 := E<ji2(0p2,(Jip2) 



E\2 := E84-cos(0p2) cos(({ip2) - E<(4 sin( 4>p2 ) 



Ey2 := E84 cos(0p2) sin(<J>p2) - E<J»4 cos(((>p2) 



Ez2 - E04 sin(0p2) 



/ex: 1 ! 

oa2 := 1 E>2 

|Ex2| ) 2 -c ( |Ey2| ) 2 -t (|Ez2|) 2 Wj 



ca2 



1 0 1 2034 + 0 48533.1 ' 
0 12034 + 0 4S533j 
r0 60556 - 0.36501) , 



PLF2 : 



2 := ( |c\v2 aalj) 



PLF2 =0.25003 



( radians ) 



( radians ) 



(V/m) 



(V/m) 



(V/m) 



(V/m) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 
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PARABOLIC REFLECTOR KITH HELICAL FEED FAR-FIELD ELEVATION PATTERNS 



E -FIELD RADIATION PATTERN : 



ll :=20 



e: 




— - io ' 6 — - io ‘ 6 

36 36 



X2( 02 ) 





1 -\ 




1 -\ 


E02 1 


H 


o 

o 

V) 





Y2(02 ) := 



E62 f 02 



sin[ 02 + — 




— E FIELD 



X2( K ? ) 
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E-FIELD RADIATION PATTERN: 



1 / 


. , ! -\ 


0 ) 1 co> (C + - 


V.'i 02 = Ec2i tC.O i sin U'2 t- - 


\ 2 / 


l 2 / 




-ooo? -ooo: -ooo: -oooi -oooi o oooi oooi ooo: 0002 000? 



— E FIELD 
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THE SrlPTT ANTENNA (FEFDECTO? OPTION 
MATH CAT' SOFTV.'ARE - DI SH_S PI .MOD 

Spiral antennas are a family of two or three dimensional devices that 
possess frequency independant parameters over a wide bandwidth. Spiral 
antennas are commonly used for direction finding, satellite tracking 
and missile guidance. 

The planar spiral may be of the Archemedean, log-spiral, or equiangular 
type. All three radiate two main, circularly polarized lobes 
perpendicular to the plane of the antenna. Additional gain for planar 
spirals may be achieved by placing a metal cavity on the side of the 
antenna with the unwanted lobe. The cavity may be empty or be filled 
with electromagnetic energy absorbing material. These applications 
pnnciply examine the equiangular planar spiral and do not account for 
cavity backed effects. 

(Note: Math cad equations cannot use symbolic subscripts. Therefore, 

symbols like /. will immediately fellow the parameter m equations, ) 

The spiral antenna Mathcad applications will compute the following 
parameters for equiangular planar spirals: 



k = Wavenumber 
/. = Wavelength 

D 0 = Directivity 

E<t’ - Electric Field Component 

U = Radiation Intensity 

Ujnax = Maximum Radiation Intensity 

p rad = Radiated Power 

G ~ Ga i n 

EIRP = Effective Isotropic Radiated Power 
A em = Maximum. Effective Aperture 
BW = Bandwidth 

r min = Minimum Distance to Fa r- Field 
R r = Radiation Resistance 
h ern = Maximum Effective Height 
f high = Upper Operating Frequency , 

^low = Dower Operating Frequency 

r n = Any point on the n th edge of a spiral 
e ex = Equiangular Planar Spiral Expansion Ratio 
Zi = Planar Spiral Input Impedance 
T = Voltage Reflection Coefficient 
e rv = Reflection Efficiency 
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~ c c 



. OW 



= Upper Opera tire V.'aveler.c::. 
= Lower Ope ratine Wavel enot r. 



A = Planar Spiral Elect: 



Fielc Amplitude 



The followmc data must be incut based on known cr estmated a a t a : 



M = Mode 

N = Nuirrer of Spiral Arms 
f = Frequency of Interest 

1 = Number of Increments for Far Field Radiation Patterns 
r^^ = Distance of Far-Field Calculations 

i. i. 

I c = Input Current at Antenna Terminals 
r G = Spiral Feed Point 
a = Flare Rate 

6 n + i = Angular Arm, Width of n tn Spiral Arm 
$ r = Azimuth to Compute Expansion Ratio 
P = Overall Radius 

E c = Source Strength Constant for Planar Spirals 

c w = Wave Unit Polarization Vector 

Z 0 = Characteristic Impedance of Feed Assembly 

o - — Equiangular F 1 a n a r Spiral Unit Polarization Vector 



THE PLANAR SPIRAL ANTENNAS 



r 



er 



data here: 



N =2 
M = 1 



( a rms 



(mtooe 



f = 3 5 10 10 (Hz. 



lo := 1 



( amps , 



{Note: K n . a> . is N-l: 
l •= 30 (increments' 



rff : 1 10' (meters 



a ‘= .221 ( dimensionless ; 

R - 01 (ml 



ro = 00 1 { m ) 



Eo := 10 ? 



(V/m) 
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’ 0 




1 


2 




ji 




( radi ans ou := 


1 


71 




7 


3 z 




y- 


2 




0 



(dimensionless ) 



Zo := 100 



(Q 



Calculate planar spiral antenna geometric parameters and define 
constant s : 



c:=2 9979 10* (meters / sec) 
c 

a =- (meters / cvcle 

f 

a = S 5654?* 10 ' (meters / cycle’ 



oa 




(dimensionless ) 



V* 120a 



(Q) 



. 1 in -9 

c 0 • T - ICr 

j> 0 n 



(Farads / m) 



■ = 4 r. 10' 



( H / m ) 



Calculate 



planar spiral antenna parameters : 



Define angular offset 9 
from v-z axis: 



e 




i 



-- io' 6 1- lcr 6 

2 2 



( radians ) 



*> 7 

<t> r 0 12-.' 2n (radians) 

l 
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Lis 



a r.:e 



r c r - : 1 : 



rr 0 '= 1 6 ). { r - 



n ] = 1 t < R 



8 R : 




rmin ■= ma\( rr) (m / 



rmin = 0 1 



Wavenumber k: 



~ i stance 



n±J2 Sn : 



*■» „ 
/. 



k = 7 33552 1 0‘ rr - : p i 2 r j = 4 ( mjv p* 10 



v n c: = ro e 



4-P 



i rr 



Expansion F.atic c ^ > . j_ 



cexf n.$r) 



v( n . er «h 2 r ! 



r( n. cr ) 



ci n. 0 r. s 



ess 



ee\( 1 ,2 a) =4 0091" 



ess 



Bandwidth BV,’ : 



Ecui ar.auia r Snral 



joe- Periodic Snral 



/.high := 4 ro 



/.high 1 := 20 ro 



(rr. 



fhigh : = - 



/.high 



(Hr 



fhigh 1 ■ = 



/.high! 



(Hr - 
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lhigh = ” 4 94 "5 10 



fhizh = " 4v4~5 10 



/.low := 4 R 



(m 



flow *= • 



/.low 



(Hz) 



How = 7 49475 10 



(H: 



BW .= fhigh - flow 



(Hz 



BWl fhigh 1 - How 



(Hz) 



BW = 6 74528 10 



10 



(HZ 



BWl =7 494^5- 10 



(HZ' 



Electric Field £( 6j : and Electric Field Amplitude A( g ;: 



w =0 i (increments) 



M 



7 M m j 

Eo k‘ cos(0) ( 1 4 j acos(G)) lan - e 



M I r + - j — k rff 



E(0) :=- 



sin(0'f rff 



( V / rr. / 



A(0) =- 



cos(0)tan|-| e 



M I—) aum a cos>( 0 j ) 



sin(0j 



/ j 

1 + a fc cost 0) 



V/rr. 



H'( 0) 



M { 2 \ 

= — ln\l + a" cos(0) J + atan{ a cos(G) ) 



(dimensionless ) 



E4>( 6) := A(0) e* J 



J 

• 4 i( 0 ) e 




rff 



(V/m) 
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(V/m) 



A1 :=- 



.71 w , 

cos| - - + 71 — | tan 



71 W 

- + 71 — 

2 i 



M 






r 7i \\ ) j . 2 r 7i \\ 

sm|---+7i — M-+a C0S1-- + 7 — 

^ i N 



Radiation Intensity U ( 6 ) : 

U(0) :=— (a(0) 2 ) 

2t 1o 




Umax := max(Ul ) 



Umax =2.2285-10" 3 



Radiated Power P rarj : 

'n 



Prad n 



U(0) sin(0) d0 

0 



Prad = 7 94528* 10”' 



(W / solid angle) 



(W / solid angle) 



(W / solid angle) 



(W / solid angle) 



(V? 



(W) 



Directivity D n : 

^ 4n Umax 

D = (dimensionless, 

0 Prad 

D 0 = 3.52462 (dimensionless) 
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Rr •= 



2 Prad 

(|lo |) 2 



(Qj 



Rr =001589 (fi; 



(Q 

, M 

sin | n — 

N 

Zi = 1 88496- 10 J (Q) 



Vci taae Reflection Coefficient F : 



r := 



Zi - Zo 
Zi + Zo 



( dimensionless 



Reflecton Efficiency c r ,,: 

cn := 1 - drlr (d imensicnless 



f = U 306"5 (dimensionless ot =0 90591 ( oimensi onl ess j 



Gam G i 



G '= cn D q (dimensionless ) 



GdB := 10 log(G) (dB) 



G = 3 1 9298 (dimensionless 1 



GdB =5 04196 



(dB 



Effective Isotropic Radiated Power t'EIRP': 

Prad G 



EIRP :=« 



(W) 



cn 



EIRP =0 028 (W) 



Polarization 




Loss 



2 



Factor ( PLF ; : 



(dimensionless ) 



PLF = 1 



(dimensionless) 



194 




a : D 0 eh PLI 

Aem ’ = (r** 

4 7 



Aem = 1 86416 icf 5 (m : 



Maximur. Effective Height : 



hem 






Rr Aem 

% 



(rr. 



hem 



5 0063* ](" J 



. r. 
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THE EQUIANGULAR PLANAR SPIRAL ANTENNA FAK- FIELD ELEVATION PATTERN 



For the purpose of this far-field radiation pattern, the spiral antenna 
lies parallel to the Ey = 0 grid line and is centered at the origin. 

The magitude of the electric field pattern is rotat 1 onally symmetric 
with respect to the Ex = 0 grid line. The equiangular planar spiral 
antenna possesses a mirror image radiaton pattern m the -y half plane. 

E<t>x( G) := | E*( 0) | E$y(0) := |e<k 6)| smlO-fL 




— E FIELD 
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PARABOLIC REFLECTOR ANALYSIS WITH SPIRAL ANTENNA FEED 



The spiral antenna Mathcad application (Reflector Option) will compute 



the following ada:t:cr,al parameters (Note: a 2 or 22 will follow most 



secondary field parameters m this application, and inorder tc speed up 



comput 


:ation time corresponding vectors are created with subscript A 6 B 


p x,y 


= E-Field Integral and Summation Equation Vector Components 



E a >:,ay = Aperture E-Field Components 



e\ P = 


= E-Field Integral and Summation Equation Parameters 


E e = 


Secondary Field Theta E- Field Component 


E * = 


Secondary Field Phi E-Field Component 


u = 

p rad 


Secondary Field Radiation Intensity 
= Secondary Field Radiated Power 


D 0 = 


Secondary Field Directivity 


G = 
c ap = 


Secondary Field Gain 
Aperture Efficiency 


c s 


Spillover Efficiency 


c t 


Taper Efficiency 


c p = 


Phase Efficiency 


c x 


Polarization Efficiency 


c ohmi c 


= Ohmic Efficiency 


A e = 


Reflector Effective Area 


^e 


Reflector Effective Height 


EIRP 


= Reflector Effective Radiated Power 


E x , y , : 


7 = Secondary Field Cartesian Components 


°a = 


Reflector Unit Polarization Vector 


PLF = 


= Reflector Polarization Loss Factor 


r min 


= Minimum. Distance to Secondary Far-Field 


CD 

o 

It 


Reflector Periphery Coaltitude 


W 


= Theta and Phi for Polarization Loss 



The foliowinc additional data .must be inputted based on known or 



estimated data 



a ref 


= Radius of the Mcutn of the Reflector 


foe = 


= Reflector Focal Length 


E b = 


Blockage Efficiency 


c sp 


Spar Efficiency 



1 1 , 1 2 - Radiated Power Increments for Phi and Theta 
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ll = Secondary Field Increments 
N, M = Summation Increments 

r ff2 = Secondary Field Observation Distance 
x, y, z = Polarization Loss Factor Coordinates 
o w = Incoming Wave Unit Polarization Vector 



Enter additional input data here: 



aref := 0.32 



(meters ) 



foe := 0 147 



(meters ) 



cb •'= 96 



esp := .95 



(dimensionless ) tl :=10 



t2 := 10 



(increments ) 



il := 10 



N := 10 



(increments ) 



M := 10 



rfD := 10 ? 



(meters ) 



x •'= 1000 



crvv2 



0 



y := 1000 



( dimensionless ) 



z = 1000 



( dimensionless ) 



(increments) 



( increments ) 



( increments ) 



(m) 



(m) 



(m) 
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Establish Integral 



arb Summa tier: Inc r erne n t s : 



*1 := 0 ,— 2 * 

tl 



radi ans ) 



A = ) N 



( increments ) 



42 •'= \(J~ t — «t 1 CT ' 2 ::+ 1 CT' (radians) 
1 ] 



B = 1 M 



(increments ) 



J fl \ ^ "I r ~ -I 

G 2 :=- — - icr,— I - icr-\ __ jo- 

136 / il 18 36 36 



♦n B *2 # — 

B M 



( radians ) 



rp := 0 aj-ef ( rr, ) 

t2 



$21 := 2 * 



. = , ^ 



M 



( radians ) 



rp ) := aref — 

A N 



(m) 



022 ■= [ — 1 — (radians) 

A \ 18 / N 36 



e: 2 o : = ’~ 



(radians ) 



Calculate reflector geometric parameters: 



. . / aref \ 

0o •= 2 aian 

\2 foe ) 


( radians ) 


— 1 
0 

11 

ON 


(rr.) 


0o = 1 65544 


(radians 


r = 10 aref 


(m) 


4 foe* + rp 2 

pCrp) 

4 foe 


(rr.) 


8 aref 

V = . 


(m) 


Op(rp) := 2 atari | — — — ] 
\2 foe/ 


( radians ) 


nmin '= max(r) 


(m) 


4 foe 2 + (rp ] A ) 2 
P '*" 4 foe 


(m) 


rmin = 95.64028 


(m) 


/rpl A \ 

Op] A ’= 2 atan 

A \2 foci 


( radians 
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Aperture Electric Fields E pv , E a v j_ 



Ea\’(rp,$l ) := E<K&p(rp)) sin(<M ) e 



I k ( j>< rp ) — fou 2 ) 



(V/m) 



Eay(rp,*l) =E4(ep(rp)) cos(4il) e* 1 ^ > - (oc 2 > 



(V/m) 



£3X^,3 sin(*ll B ) e ‘'^ lA-r ° C ^ 



(V/m) 



E ayl A B : = E *(epl A ) cos(*ll B )e 



.1 •l' (plA - foc 2 ) 



(V/m) 



Electric Field Integral Equation Vector Components P x> y j_ 
’ aref 



Ax(62,42,*l ) •= 



Eax(rp,*]) e' k n ’ sin<02) (cos(C1) 42 >>rp ^ 



0 



(V/m) 



Px(02,4>2) •= 



Ax( 02 , 4*2 , 4> 1 ) d$l 



(V/m) 



Ay( 02 , 4>2 , 4> 1 ) •= 



aref 



Eay( rp , <J> 1 )-e > ^ n> ^,92, icostti) cos ( «) + s,n<«ii sm<c2)) 



(V/m) 



Py(02,*2) •= 



2 71 



Ay(62,*2,*l )d*l 



(V/m) 



Axl 



Eaxl c krpls,n ( e22 - J ) ( C0S ( tn B) cos(e22 E )- l . S in(eil B ) $ . n («2 B )) (V/m) 

Jo 



Pxl 



A.B 



1 7! 



^A^ 11 



(V/m) 



Avl 



A,B 



Py] 



A,B 



aref 

0 

2 71 



Eay\ AB e krpls,n ( 022 ^) ( cos ( 01, b) c< *(« 22 b) sm(o22 B )) ^ ^ (V/jn) 






(V/m) 



0 



2 00 



(V/m) 



Secondary Field Electric Field Components E ft, 



E02(02,*2) =i-Ei 



,-J k rfH! 



2 7t rfi2 



( P\( 62 , *2 )• cos( i2 ) + Py( 02 , 42 ) sin( *2 ) ) 



E$ 2 (e 2 ,* 2 ) = 



. J -k-e’ 



j k rir: 



2 7i rfi2 



(P\( 02,^2) - sin($2 )+ Py(02,$2) cos( $2 ) ) cos( 02 ) 



* -j k rff? 

E63. „ =i— : (PxI 4 b c «(*::b)* p >1a.b si "(« : b)) 



A ’ B 2 7i rfT2 



, -j k rfT2 

« 3 ab r.rfn h'».B “"(*“>)) * ^‘a.b “ s (*-b)] c “( k -J 



Secondary Field Radiation Intensity U 
i .' 2 



U2(62,i2) :=— [( |Ee2( 02 .42)1 ) 2 + ( | E42f 02 .*2 )| ) 2 ] 

‘ ^ 



( o 



U3 ab : ~[(I e 93 a.bI) : *(I e «abI) : 

2 T ) 0 



Secondary Field Radiated Power P rar j: 

Prad: * (tt) (ttj) -K 



Prad2 = 1.84163 10 



“ii 



Secondary Field Directivity D n : 

D o 22 := 4 BmaxrU?) 

Prad2 



Do22 = 4.59318 10 5 



(V/m) 



(V/m) 



(V/m) 



(W/solid ang) 



(W/s olid ang) 



(W; 



(W) 



(dimensionless) 



(dimensionless ) 
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( dimensionless ) 



Spillover Efficiency e c : 



es 



Go 



ict 



U(Gp) sin(0p) dOp 



icr 6 



U(Gp) sin(Gp) dOp 






es = 1.00042 



Taper Efficiency : 



fGo 






2 coil [£ 



Go 



HMU' slM 8p)dep 

Prad 



et =045135 



Phase Efficiency e^ 

» 4 t 6 10 " a ref 

X 

cp •= e 

Ep =0.97222 



Aperture Efficiency c a p j_ 

ex *- .98 (dimensionless ) eohmic •- 98 



eap .= es d ep-ex eb -esp Eohmic (dimensionless ) eap = 0.38452 



(dimensionless ) 



(dimensionless ) 



( dimensionless ) 



(dimensionless ) 



( dimensionless ) 



(dimensionless ) 



(dimensionless ) 
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a tea 



Effective Isotrcp;c Radi 
EIRP22 := Prad2 Do22 

EIRP22 =8 45893-1 O' 6 

Gain G: 

G22 '= cap Do22 (dimensionless 

G22 = 1 76615- 10 s (dE 

Effective Aperture : 

*> 

Ae22 *= cap ti arer 
Ae22 =0.1237 

Radiation Resistance P .: 

Rr22: S 2-^ r . 

(I iolr 

Rr22 = 3 68326 l(f u 



(VD 



(W) 



G22dB = 10 log( G22 > ( dimensionless j 



G22dB = 52 47028 (dB) 



(m 2 ) 



(m^ ) 



(Q) 



(Q) 



Effective Heiqht h f 



he22 •= 2 Rr22 



Ae22 



he22 =2 19868 10 



r 7 



(m) 



(m) 
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Polarization Loss Factor PLF: 



ep2 




( radians ) 




6p2 =0.95532 



( radians ) 



4>p2 = 0.7854 



E04 :=E62(ep2,(t>p2) (V/m) 



E<H := E()>2(ep2, 



Ex2 E04 cos(0p2) cos($p2)- E<j4 sin(<Jp2) 



Ey2 E94 co$(0p2) sin(4>p2) - Ef4 cos($p2) 



Ez2 := -EB4 sin(0p2) 



ca2 = 



1 

J(|Ex2|) 2 4(|Ey2| ) 2 4 (|Ez 2 |) 2 



/Ex2\ 

E\2 

\bz2I 



L 0 .40672 4 0 03531.) 



oar = 



-0 40672 + 0.0353 lj 



0 81 344 - 0.0706 3j 



PLF2 : 




aw2 




PLF2 =0.16667 



( radians ) 



( radians ) 



>2) (V/m) 



(V/m) 



(V/m) 



(V/m) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 
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PARABOLIC REFLECTOR WITH SPIRAL FEED ANTENNA 



E-PANE RADIATION PATTERN: 
il := 20 



X 1 ( 62 ) 



E62 1 62 



62 .= - 




cos 




10 ’\- 



_L_ icT. — - 10' ? 



IS 36 



36 



Y 1 ( 62 ) '= 



E62|62,2 



sm| 62 -+ — 




— E FIELD 
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H-PLANE RADIATION PATTERN: 



Y2/02) 



X2(02) := |e*2(02.O)| 



cos| 02 + — 



Y2(02) := |e^ 2( 02 ,0)| sin^02 




— E FIELD 



X2( 02) 
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THE PYRAMIDAL HORN ANTENNA (REFLECTOR OPTION) 

MATHCAD SOFTWARE-DISH_PYR.MCD 

Pyramidal horn antennas are devices used to provide a transistion from a 
rectangular waveguide to an unbounded medium such that the wavefront at 
the aperture of the horn has nearly a constant phase at any point in the 
mouth of the horn. 



Pyramidal horns are the most popular type of feed elements for 
reflectors used in satellite tracking, microwave communications, and 
radar . 



The pyramidal horn applications may be used to analyze E- and 
sectoral horns. To analyze an E-plane sectoral horn set horn 
(al) equal to waveguide dimension (a) . To analyze an H-plane 
horn set horn dimension (bl) equal to waveguide dimension (b) 



H-plane 

dimension 

sectoral 



(Note: Mathcad equations cannot use symbolic subscripts. Therefore, 
symbols like l will immediately follow the parameter in equations in 
lieu of subscripts.) 



The pyramidal horn antenna Mathcad applications will compute the 
following parameters: 



k = Wavenumber 
/. = Wavelength 

D Q = Directivity 

Eq = Electric Field (0) Component 

E^ = Electric Field (<()) Component 

U = Radiation Intensity 

Umax = Maximum Radiation Intensity 

P ra d = Radiated Power 

G = Gam 

EIRP = Effective Isotropic Radiated Power 
A em = Maximum Effective Aperture 

BW = Bandwidth 

r min = Minimum Distance to Far-Field 

R r = Radiation Resistance 

hem = Maximum Effective Height 

E x , y f z = Electric Field Components m Cartesian Coordinates 
°a ~ Unit Polarization Vector 

f c = Transverse Electric and Transverse Magnetic Cutoff Frequencies 

G ap = Aperture Efficiency 

PLF = Polarization Loss Factor 
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Dp = Pyramidal Horn Corner to Corner Distance 

p e , h = Pyramidal Horn Perpendicular Flare to Mouth Distances 
Il,2 = Electric Field Component Equation Coefficients 



The following data must be input based on known or estimated data: 



t = Number of Cutoff Frequencies Calculated 
m, n = Modes 

f = Frequency of Interest 

i = Number of Increments for Far Field Radiation Patterns 
rff = Distance of Far-Field Calculations 
I Q = Input Current at Antenna Terminals 

0g = Coaltitude (Deflection Angle from +z Axis) for Azimuth Plot 
<j)g = Azimuth Angle for Elevation Plot 

(x,y,z) = Coordinates for Unit Polarization Vector 

a,b = Rectangular Waveguide Dimensions 
a l'bi “ Pyramidal Horn Dimensions 

Pl'P2'Pe'Ph = Pyrandcial Horn Imaginary Cone Apex to Mouth Distances 
E q = Electric Field Amplitude Constant 



Enter input data here: 



t -= 5 (modes ) 



f ■= 3 1 0 8 (Hz) 



i - = 1 8 (increments) lo : = 1 (amps) 

rff — 110" (meters ) 



(Note: For E-plane sectoral horn analysis 
set (al) equal to (a) . For H-plane sectoral 
horn analysis set (bl) equal to (a)) 

a := .5 (m) b := .25 (m) 

al .= 12 (m) bl •= b (m) 
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pi — 6 



(m) 



<nv 



±_ 

6 



(dimensionless ) 



p2 := 6 
\ •= 10 3 



pe 6.15 (m) 

ph:=6.6 (m) 



y := 10 3 
z — 10 3 



Eo := 1 (V/m) 



Calculate pyramidal horn antenna geometric parameter. 
constants : 



c -= 2.9979- 10 8 (meters/sec) 



TU— 120-ti 



c 

K •=- (meters/cycle) 

f 



- 1 WT 9 

c ict 

0 36 a 



X = 



(meters/cycle ) 



u 0 : = 4 n 10" 



Dp 



:= + 



bl (m) 



Dp = 



(m) 



pe := (bl - b) 




(m) 



ph := (al - a) 




pe = 



(m) 



ph = 



(m) 



(m) 



(m) 



(m) 



and define 



(Q) 



( Fa rads/m) 



(H/m) 



(m) 

4 



(m) 
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Calculate planar spiral antenna parameters : 



Define angular offset Q 
from y-z axis: 



Distance to Far-Field r m in : 



0 : = 



7171 71 , 6 71 1A .6 

+ 1(T + 10 

2 i 2 2 



( radians ) 



IT =16 1 



n'j := 5 Dp 



4> := ia 6 ,— + 10' 6 .. 2 jh- 10’ 6 
i 



( radians ) 



IT_ 




(m) 

(m) 



(m) 



rmin := ma\(rr) (m) 

nnin = (m) 



Wavenumber k: 



k =— (m 1 ) 

1 

k= (rrT 1 ) 



Cutoff Frequencies f ^ : 



n := 0.. t 



(modes ) 



m := 0 t 



(modes ) 




(Hz) 
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(Hz) 



fc = 



(Note: The index for both m, n above begins with zero, 
have m or n equal zero. TE qq mode does not physically 



Electric Field Component Coefficients 
v := 0 . i 

\v = 0. 1 



11(0.4*) 



' al 
1 



al 

i 



cos 



M) 

\ V 



l- 2 

a -J k — - 

\2 pi 



sin( 0 ) oos( o ) ' 



' d! 




11(6,4) =if 



1- a 


cos|k - sin(0) cos( <fr)J 


,11(6,4) 


2 


(k — sin(0)cos(4) 


2 -[A 7 




[\ 2 J 


\2j \ 





TM modes cannot 
exist . ) 



(increments ) 
(increments ) 



(dimensionless ) 



(dimensionless ) 



( dimensionless ) 
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